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ABSTRACT 

We analyze the ASCA spectra accumulated within ~ 100 kpc radii of 12 of the 
brightest groups of galaxies. Upon fitting isothermal models (IT) jointly to the ASCA 
SIS and GIS spectra we obtain fits for most groups that are of poor or at best marginal 
quality and give very sub-solar metallicities similar to previous studies, (Z) — 0.29 ± 
0.12 Zq. Two-temperature models (2T) provide significantly better fits for 11 out of 
the 12 groups and in every case have metallicities that are substantially larger than 
obtained for the IT models, (Z) = 0.75±0.24 Zq. Although not very well constrained, 
for most of the groups absorption in excess of the Galactic value is indicated for 
the cooler temperature component of the 2T models. A simple multiphase cooling 
flow model gives results analogous to the 2T models including large metallicities, 
(Z) = 0.65 ± 0.17 Zq. The nearly solar Fe abundances and also solar a/Fe ratios 
indicated by the 2T and cooling flow models are consistent with models of the chemical 
enrichment of ellipticals, groups, and clusters which assume ratios of Type la to Type 
II supernova and an IMF similar to those of the Milky Way. 

Thus, we have shown that the very sub-solar Fe abundances and Si/Fe enhance- 
ments obtained from most previous studies within r ~ 100 kpc of galaxy groups are 
an artifact of their fitting isothermal models to the X-ray spectra which also has been 
recently demonstrated for the brightest elliptical galaxies. Owing to the importance 
of these results for interpreting X-ray spectra, in an appendix we use simulated ASCA 
observations to examine in detail the "Fe bias" and "Si bias" associated with the 
spectral fitting of ellipticals, groups, and clusters of galaxies. 

Key words: galaxies: general - galaxies: evolution - X-rays: galaxies. 



1 INTRODUCTION 

The hot gas in elliptical galaxies, groups, and clusters of 
galaxies is thought to be enriched with metals by material 
ejected by supernovae. Hence, the abundances in the hot gas 
record the star formation history of these systems and there- 
fore probe the rates of Type la and Type II supernovae as 
well as the shape of the IMF (for reviews see Renzini 1997, 
1999). Since the hot gas in these systems emits at tempera- 
tures T ~ 10 7 -10 8 °K (~ 1-10 keV), X-ray observations are 
the preferred means to determine the abundances in the hot 



Spatial analyses of the X-ray s urface brightnes s profiles 
of grou ps of galaxies using ROSAT ( Triimper 198ot) indicate 



that th i hot gas typically extends out to distances in excess 
of 100 kpc (e.g. David et al. 1994; Pildis, Bregman, & Evrard 
1995; Davis et al. 1996; Ponman et al. 1996; Trinchieri, Fab- 



biano, & Kim 1997) and also suggest that the gas is divided 
into distinct components corresponding to the halo of the 
central galaxy and an extended component for the surround- 
ing group (Mulchaey & Zabludoff 1998; also see Ikebe et al. 
1996). These studies typically obtain average temperatures 
of T ~ 1-2 keV from analysis of the ROSAT PSPC spectra 
and also find evidence for positive temperature gradients in 
some groups (Trinchieri et al. 1997; Mulchaey & Zabludoff 
1998). Despite the evidence for temperature gradients, in 
most cases the abundances reported by the ROSAT studies 
are obtained from isothermal spectral models from which 
very sub-solar Fe abundances are inferred (a notable excep- 
tion being NGC 533 - Trinchieri et al. 1997). 



Similarly, analyses of the ASCA data (Tanaka et al 



1994) of groups which assume isothermal gas typically find 
very sub-solar Fe abundances and, in some cases, Si/Fe ra- 
tios in excess of solar (Fukazawa et al. 1996, 1998; Davis, 
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Mulchaey, & Mushotzky 1999). Although the point spread 
function (PSF) of ASCA is quite large and asymmetric, 
there have been several attempts to deconvolve the ASCA 
data to allow analysis of the spectra as a function of ra- 
dius with the results for galaxy clusters being decidedly 
mixed (e.g. Markevitch et al. 19 98; Kikuchi et al. 1999). 
Recen tly, F inoguenov et al. ( |l99S| ) and Finoguenov & Pon- 
man (199!:) have applied their own method of deconvolu- 



tion to the ASCA data of a small number of groups. By 
assuming the hot gas to be in a single phase, these authors 
have fitted isothermal models at several radii for each group . 
Similar to the previous stu dies, Finoguenov et al. (199£) 
and Finoguenov & Ponman obtain very sub-solar Fe 

abundances and Si/Fe ratios in excess of solar. 

The very sub-solar Fe abundances for the hot gas in 
groups of galaxies deduced by these studies are in many cases 
less than those of galaxy clusters and thus pose a problem 
for the popular models of chemical enrichment which as- 
sume ratios of Type la to Type II supernova and an IMF 
similar to those of the Milky Way. (e.g. Renzini 1997). To 
reconcile these observations with these models it has to be 
postulated that galaxy groups accrete a significant amount 
of primordial gas after they have spent most of their evolu- 
tion expelling gas. The un attrac tive features of this scenario 
are discussed by Renzini (1994). 



These results for groups are in every way analogous 
to those initially obtained for ellipticals with ROSAT (e.g. 
Davis & White 1996) and ASCA (e.g. Matsumoto et al. 
1997). Recently, however, Buote & Fabian (1998, hereafter 
BF) and Buote (1999, hereafter B99) have demonstrated 
that models having at least two temperature components 
in the hot gas are required to explain the ASCA spectra of 
the brightest ellipticals. These multiphase models agree with 
the positive temperature gradients indicated by the ROSAT 
data unlike the isothermal models (B99). 

Perhaps most importantly BF and B99 also showed that 
the Fe abundances determined from the multiphase spectral 
models of the brightest ellipticals are consistent with so- 
lar in stark contrast to the very sub-solar values obtained 
from isothermal models. This strong dependence of the Fe 
abundance on the spectral model indicates that one cannot 
obtain a reasonable estimate of the Fe abundance by us- 
ing an isothermal model if the hot gas consists of multiple 
temperature components. Similar problems also arise when 
determining the a/Fe ratios (B99). 

Since the spectral properties of the brightest ellipticals 
are so similar to groups, it should be expected that similar- 
results for the Fe abundance also apply to groups of galaxies. 
To investigate this issue we selected for an alysis the groups 
from t he sam ples of Mulchaey & Zabludoff (1998) and Davis 
et al. ( 1999| ) having the highest quality ASCA data suitable 
for detailed spectral analysis. In addition, we inc luded the 
group NGC 1132 since Mulchaey & Zabludoff (Il999|) ar- 



gue that it has properties similar to bright groups despite 
the fact that it has no large galaxies aside from NGC 1132. 
Our final sample consists of 12 galaxies whose properties are 
listed in Table |. 

The paper is organized as follows. In section [| wc 
present the observations and discuss the data reduction. We 



define the models used for broad-band spectral fitting^ in 
section |3.1.l| and discuss calibration and fitting issues in 
section p. 1 .2) . The results for the broad-band spectral fitting 
analysis are presented in section We discuss the implica- 
tions of our results in section ^| and present our conclusions 
in section^. Finally, in Appendix ^|we use simulated ASCA 
observations to examine in detail the "Fe bias" and "Si bias" 
associated with the spectral fitting of ellipticals, groups, and 
clusters of galaxies. 



2 ASCA OBSERVATIONS AND DATA 
REDUCTION 

The ASCA satellite consists of two X-ray CCD cameras 
(Solid State Imaging Spectrometers - SISO and SIS1) and 
two proportional counters (Gas Imaging Spectrometers - 
GIS2 and GIS3) each of which is illuminated by its own 
X-ray telescope (XRT) . The SIS has superior energy resolu- 
tion and effective area below ~ 7 keV while the GIS has a 
larger field of view. Although the PSF of each XRT has a 
relatively sharp core, the wings of the PSF are quite broad 
(half power diameter ~ 3') and increase markedly for ener- 
gies above a few keV (e.g. Kunieda et al. 1995). As in our 
previous studies (e.g. B99) we do not attempt to analyze 
the spatial distribution of the ASCA data of these sources 
because of the large, asymmetric, energy-dependent PSF. 
Rather, we analyze the ASCA X-ray emission within a sin- 
gle large aperture for each group which encloses the region 
of the most significant emission (see below). 

We obtained ASCA data for the 12 groups from the 
public data archive maintained by the High Energy As- 
trophysics Science Archive Research Center (HEASARC). 
The properties of the observations are listed in Tables ti 
and pi We reduced these data with the standard FTOOLS 
(v4.2) software according to the procedures described in The 
ASCA Data Reduction Guide and the WWW pages of the 
ASCA Guest Observer Facility (GOF)Q 

2.1 SIS 

Most of the SIS data for these groups required corrections 
and screening in addition to those performed under the stan- 
dard Revision 2 Data Processing (REV2) . Many of these ob- 
servations were taken after 1994 at which time the effects of 
radiation damage on the SIS CCDs began to become impor- 
tant (see Dotani et al. 1995). One of the problems caused by 
radiation damage is the Residual Dark Distribution (RDD) 
which is es sentially an increase in the dark current of th e 
CCDs (see jittp: / /www. astro. isas.ac.jp / ~dotani /rdd.htmj ) . 
RDD degrades the energy resolution and the detection effi- 
ciency of the data. 

Although in principle one can correct for the RDD ef- 
fect, the presently available RDD maps and software (i.e. 



We focus on the broad-band spectral analysis since B99 found 
that individual line ratios of Si and S place only weak constraints 
on the spectral models of the brightest elliptical galaxies which 
have data of significantly higher S/N than the groups in our 
sample. 



t See [http:/ /lioasarc. gsfc.iiasa.gov/docs_/asca/abc/abc. html and 
littp: / /heasarc. gsfc.nasa.gov/docs/asca^ 
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Group and BGG Propert 


















Group 


z 


N gal 


a 


N H 


logio L x 


BGG 


Type 




logio L B 








(km s" 1 ) 


(10 21 cm~ 2 ) 


(erg s" 1 ) 


Name 






(erg s" 1 ) 


HCG 51 


0.0258 




240 


0.12 


42.40 


NGC 3651 


El 


14.49 


43.18 


HCG 62 


0.0146 


15 


376 


0.27 


42.20 


NGC 4761 


SO 


13.30 


43.16 


MKW 9 


0.0397 




336 


0.39 


42.73 


UGC 09886 


SO 


14.29 


43.63 


NGC 533 


0.0181 


36 


464 


0.32 


42.31 


same 


E3 


12.22 


43.78 


NGC 1132 


0.0232 




253 


0.48 


42.48 


same 


E 


13.01 


43.68 


NGC 2563 


0.0163 


29 


336 


0.43 


41.76 


same 


SO 


13.01 


43.37 


NGC 4104 


0.0283 




546 


0.17 


42.52 


same 


SO 


12.93 


43.88 


NGC 4325 


0.0252 


18 


265 


0.23 


42.91 


same 


El 


14.41 


43.19 


NGC 5129 


0.0233 


33 


294 


0.18 


42.09 


same 


E 


12.89 


43.73 


NGC 5846 


0.0061 


20 


368 


0.42 


41.90 


same 


E0-1 


11.13 


43.47 


NGC 6329 


0.0276 






0.21 


42.07 


same 


E 


13.62 


43.59 


RGH 80 


0.0370 




448 


0.10 


42.74 











Properties for the groups and their brightest group galaxy (BGG); i.e. morphological types, blue magnitudes, and blue luminosities 
refer to the BGG. Redshifts are ta ken fr om NED and the RC3. The number of group galaxies (N ga j) and velocity disper sions (er) are 
taken from Zabludoff & Mulc haev ( |l998| ) for six galaxies. Velocity dispersions for oth er gal axies are from Hickson et al. jl992|) for HCG 
f land 62, from Beers et al. ( 1995 ) for MKW9 and NGC 4104, from Tonry & Davis ( 198l| ) for N GC 1132, and from Ramella et al. 
( 1995| ) for RGH 80. Galactic Hydrogen column densities (JV H ) are from Dickey & Lockman ( 1990 ) using the HEA SARC w3nh tool. The 
X-ray luminosities, L x , are computed in the 0.5-10 keV energy band from the (unabsorbed) 2T models in section 3.2.2 of this paper. 
The (luminosity) distances used for these calculations are obtained from the quoted redshifts w ith H p = 70 km s _1 Mpc _1 and 
Qo = 0.3 except for the nearby galaxy NGC 5846 for which use the distance from Faber et al. ( L98£ ) rescaled for Ho = 70 km s" 1 Mpc 
— . Morphological types and total apparent blue magnitudes (-BJJ,) for the BGGs are from the RC3. The blue-band luminosities of the 
BGGs are computed from B^, using the distances as above and LgQ = 4.97 X 10 32 erg s — 1 . 



Table 2. ASCA Observations 



Group 


Sequence 


Date 


SIS 


SIS 


RDD 


Default 


SIS Temp 


SIS BKG 


Radius 


Radius 






Mo/yr 


CCD 


Data 


Correct 


Screen 


SO 


SI 




(arcmin) 


(k P c) 








Mode 


Mode 






(°C) 


(°C) 




SIS0 


GIS 


SIS0 


GIS 


HCG 51 


82028000 


6/94 


1+2 


B 


N 


Y 






T 


3.9 


4.5 


127 


146 


HCG 62 


81012000 


1/94 


2 


B 


N 


Y 






T 


3.0 


3.5 


55 


64 


MKW 9 


83009000 


2/95 


2 


B 


Y 


N 


-61.0 


-60.8 


L 


4.0 


4.4 


201 


221 


NGC 533 


62009000 


8/94 


1 


B 


N 


N 


-59.5 


-59.0 


L 


3.5 


3.5 


80 


80 




62009010 


1/96 


2 


B 


Y 


N 


-59.5 


-59.0 


L 


3.5 


3.5 


80 


80 


NGC 1132 


65021000 


8/97 


1 


B 


Y 


Y 






T 


3.8 


4.5 


111 


131 


NGC 2563 


63008000 


10/95 


1 


F 


Y 


Y 






T 


3.0 


3.4 


61 


70 


NGC 4104 


83008000 


12/94 


2 


B 


N 


N 


-60.5 


-60.0 


T 


3.0 


4.0 


107 


143 




84017000 


5/96 


1 


F 


Y 


N 




-61.0 


T 


3.0 


3.4 


107 


121 


NGC 4325 


85066000 


1/97 


2 


F 


Y 


Y 






T 


3.0 


3.4 


95 


108 


NGC 5129 


84048000 


7/96 


2 


B 


Y 


Y 






T 


3.0 


3.5 


88 


103 


NGC 5846 


61012000 


2/94 


4 


B 


N 


N 






L 


4.1 


4.5 


55 


60 


NGC 6329 


84047000 


4/96 


2 


B 


Y 


Y 






T 


3.0 


3.4 


104 


118 


RGH 80 


83012000 


6/95 


2 


B 


Y 


Y 






T 


3.6 


3.8 


169 


178 


See text in section ^ for ex 


planation of the terminoloj 


y used in this table. Additional information for some individual 


groups is 


given 



the end of the section. Distances computed as described in notes to Table 



FTOOL CORRECTRDD) do not fully restore the SIS data for 
energies below ~ 0.8 keV. The primary result is that the SIS 
data overestimate the column density by ~ 2 x 10 20 cm -2 , 
an effect which is more serious for the SIS1 (see the ASCA 
GOF WWW pages). This excess column is of similar mag- 
nitude to the Galactic columns of the groups in our sample 
(see Table §). 

The most accurate correction for the RDD effect is pos- 
sible for SIS data taken in FAINT mode. After applying 
the RDD correction the FAINT mode data is converted to 
BRIGHT2 mode data for which spectral responses can be 
generated using the standard software. Unfortunately, for 
several of the observations the exposures were approximately 
evenly split between FAINT mode and BRIGHT mode. 
Since the spectral responses of BRIGHT and BRIGHT2 
mode data are different, the data from these different modes 



cannot be combined. However, the FAINT mode data in- 
stead can be converted to BRIGHT mode data (as is done 
in the standard processing), and this BRIGHT mode data 
can be corrected for the RDD effect. 

Since the benefit from larger S/N outweighs the small 
improvement in the RDD correction for FAINT mode data, 
for groups that require the RDD correction but have a signif- 
icant fraction of their exposure taken in BRIGHT mode, we 
performed the RDD correction on the BRIGHT mode data: 
i.e. both on the real BRIGHT data and the BRIGHT data 
that is converted from FAINT mode. For other groups that 
require the RDD correction, we applied the correction on the 
FAINT mode data which is then converted to BRIGHT2 for 
further analysis. Finally, we analyzed the BRIGHT mode 
data for all groups that do not require RDD correction. (As 
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Table 3. ASCA Exposures and Count Rates 



Name 


Sequence ^ 


Exposure 


Count Rate 


Exposure 


Count Rate 










(10.-2 c £ 


q — 1^ 
S ) 


no 3 






ct S _ ) 








DID 1 


OlOU 


DID 1 


VjlOZ 


vjjIOO 




VjjIOO 


11 L V_i 01 


oZUZoUUU 


79 n 

I z.u 


71 fi 


^ 7 


1 . u 


79 n 


79 n 


2 2 


9 f{ 




m onnn 
olulzUUU 




lO 7 
lo. I 


in o 
1U. z 




OO.o 


oo.o 


. 


-1 7 


IVlxV vv y 


oouuyuuu 


9Q n 


1 £ 8 


o. ± 


o. a 


Afi 7 

4:1), / 


A(\ 7 

4:1). 1 




'L.O 


IN \jKj Doo 


so no on no 


o . a 


1 

y . i 


.'1 


0.1 


in 
iy .d 


1 Q (\ 

iy.D 


9 7 


o. u 




ozuuyuiu 


14. y 


i /i n 
14. y 


7 i 
1 . 1 


o.z 


lo.z 


1 Q O 

lo.z 


■3 9 


O. 1 


at rip 1 1 o o 

INCjO lloz 


boUzlUUU 


30.5 


30.5 


8.2 


6.1 


29.5 


29.4 


2.9 


3.6 


NGC 2563 


63008000 


45.8 


45.8 


2.7 


2.0 


51.1 


50.4 


1.1 


1.3 


NGC 4104 


83008000 


12.2 


11.6 


5.8 


4.3 


45.8 


45.8 


3.4 


4.2 




84017000 


21.8 


20.9 


5.4 


4.6 


24.7 


24.5 


2.9 


3.5 


NGC 4325 


85066000 


19.7 


19.7 


18.1 


12.7 


29.0 


29.0 


5.4 


6.1 


NGC 5129 


84048000 


31.2 


29.5 


3.0 


2.4 


34.1 


34.1 


1.1 


1.4 


NGC 5846 


61012000 


24.8 


11.5 


17.7 


9.9 


40.4 


40.4 


4.8 


5.4 


NGC 6329 


84047000 


30.9 


32.1 


3.2 


2.7 




38.3 




1.8 


RGH 80 


83012000 


39.2 


37.6 


6.0 


1.6 


45.2 


44.2 


1.9 


2.1 



The count rates are given for energies 0.55-9 keV for SIS data before 1996, 0.65-9 keV after 1996, and 0.8-9 
keV for the CIS. The count rates are background subtracted within the particular aperture (see text in 
section Q) . 



is standard only data taken in medium and high bit rates 
were used for the SIS.) 

Most of the events files were filtered using the default 
REV2 screening criteria. The most frequent departure from 
the default screening occurred because of temperature fluc- 
tuations in the SIS. The data quality, most importantly the 
energy resolution, of the SIS CCDs significantly degrades 
when the temperatures exceed ~ — 60°C, though optimal 
performance is achieved for temperatures below ~ — 61°C. 
In addition, the maps required for the RDD correction are 
defined according to 1-degree intervals in the SIS tempera- 
ture; e.g. -61°C to -62°C. The current software allows the 
RDD correction to be made with only one RDD map for 
each SIS. We defined maximum temperatures considering 
these issues. In Table [| we list the maximum temperatures 
for each SIS for those observations that required tempera- 
ture filtering. 

We also improved the data quality for a few groups 
by using more restrictive angles from the bright Earth 
(BR_EARTH) and smaller values for the Radiation Belt 
Monitor (RBM); these occurrences are listed at the end of 
this section. Our final screening of the data involved visual 
examination of the light curves for each observation and re- 
moving intervals of anomalously high (or low) count rate. 
Note only data with the standard event grades (0234) were 
used. 

The gain of the SIS varies from chip-to-chip, and the 
variation is a function of time owing to radiation dam- 
age. We have corrected the SIS data for these effects us- 
ing the most up-to-date calibration files as of this writing 
(sisph2pi_110397.fits). 

The final processed events were then extracted from a 
region centered on the emission peak for each detector of 
each sequence. We selected a particular extraction region 
using the following general guidelines. Our primary concern 
is to select a region that encloses most of the X-ray emis- 
sion yet is symmetrically distributed about the origin of the 
region; as a result, we used circles for most of our extrac- 
tion apertures. We limited the size of the aperture to ensure 
that the entire aperture fit on the detector, an issue more 



important for the SIS because of its smaller field-of-view (20' 
square for SIS vs 50' diameter for CIS). Moreover, for the 
SIS0 and SIS1 we tried to limit the apertures to as small a 
number of chips as possible to reduce the effects of residual 
errors in chip-to-chip calibration. 

We extracted the events using regions defined in de- 
tector coordinates as is recommended in the ASCA ABC 
GUIDE because the spectral response depends on the loca- 
tion in the detector not the position on the sky. The radii for 
the SIS extraction regions are listed in Table H (For NGC 
1132 we used a square box for the SIS0 data and a rectan- 
gular box for the SIS1 data. The radius listed for this group 
in Table ^ corresponds to one-half the length of the sides of 
the SIS0 box. The SIS1 region is obtained by squashing the 
SIS0 box by 10% along one of the coordinate axes.) 

We emphasize that our SIS regions cover the area of 
significant S/N and are negligibly contaminated by any ob- 
vious sources in the field. Since the X-ray emission of these 
groups is quite extended (e.g. Mulchaey & Zabludoff 1998), 
our apertures typically enclose ~ 50 per cent of the total 
emission. However, because the surface brightness declines 
rapidly with radius, and our apertures enclose the regions of 
significant S/N, the back-scattering of the emission at larger 
radii i nto o ur apertures is neglible. In particular Takahashi 
et al. (1995) showed that approximately 94% of the GIS flux 
of a model cluster within an (on-axis) aperture r = 4' ac- 
tually originates from within that aperture; an even larger 
percentage applies for SIS data because the SIS+XRT PSF 
is smaller than the GIS+XRT PSF. Since the ASCA XRT is 
smallest for energies ~ 0.5-3 keV which are most important 
for groups, the back-scattering of emission from larger radii 
into our apertures is at most a few percent and thus has a 
negligible impact on our spectral analysis. 

Our apertures differ from those of Davis et al. ( 199£ ) 



who instead defined most of their radii to encircle all of the 
emission detected with ROSAT. Their procedure leads to 
radii that are factors of 5-10 larger than our values listed 
in Table |i[ However, the ASCA data at such large radii 
are very noisy and do not improve constraints on spectral 
models obtained from using our smaller regions; i.e. we find 
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that using larger apertures increases the noise and generally 
loosens constraints on the models obtained using data in the 
apertures enclosing the highest S/N data. 

For most of the observations we computed a background 
spectrum for each detector using the standard deep observa- 
tions of blank fields; these are indicated by "T" under "SIS 
BKG" in Table |^. There are important advantages to using 
these background templates instead of a local background 
estimate. First, the templates allow background to be ex- 
tracted from the same parts of the detector as the source 
and thus the vignetting and other exposure effects are the 
same for each; this is not the case for background taken from 
a different region of the detector as the current software does 
not allow the required corrections to be made for spectral 
analysis. Second, as stated above, it is known from ROSAT 
observations of these galaxies that their emission extends 
over much of the field of view of the SIS and perhaps more 
(e.g. Mulchaey & Zabludoff 1998). 

However, it should be remembered that the background 
templates are most appropriate for SIS data taken in 4-ccd 
mode early in the mission. For more recent data taken in 
1-cdd or 2-ccd mode, using these background templates can 
lead to spurious effects particularly for energies above ~ 7 
keV where instrumental background dominates the cosmic 
X-ray background. Systematic errors in the SIS at these high 
energies are not so important for our study since the con- 
straints on spectral models are dominated by data at lower 
energies. The integrated spectra of these sources within our 
apertures are not dominated by background and, as a result, 
the deduced temperatures and abundances are not overly af- 
fected whether the templates or local background are used. 

For those sources where the screening differed signifi- 
cantly from the REV2 standard we used a local background 
estimate. The background regions were chosen to be as far 
away as possible from the group center and any other sources 
in the field. The observations using a local background are 
listed with a "L" under "SIS BKG" in Table |. 

The final screened background-subtracted exposures 
and count rates are listed in Table ^| 

The instrument response matrix required for spectral 
analysis of ASCA data is the product of a spectral Redistri- 
bution Matrix File (RMF) and an Auxiliary Response File 
(ARF). The RMF specifies the channel probability distribu- 
tion for a photon (i.e. energy resolution information) while 
the ARF contains the information on the effective area. An 
RMF needs to be generated specifically for each SIS of each 
observation because, among other reasons, each chip of each 
SIS requires its own RMF and the spectral resolution of the 
SIS is degrading with time. We generated the responses for 
each SIS using the FTOOL SISRMG selecting for the stan- 
dard event grades (0234). Using the response matrix and 
spectral PI (Pulse Invariant) file we constructed an ARF 
file with the FTOOL ascaarf. 

The source apertures used for the SIS overlapped more 
than one chip for NGC 5846 and NGC 6329. In order to 
analyze the spectra of such regions we followed the stan- 
dard procedure of creating a new response matrix that is 
the average of the individual response matrices of each chip 
weighted by the number of source counts of each chip (i.e. 
within the source aperture). (Actually, the current software 
only allows the RMFs to be averaged and then an ARF is 
generated using the averaged RMF. This is considered to 



be a good approximation for most sources.) Unfortunately, 
some energy resolution is lost as a result of this averaging. 
For the observations of our sources, however, this small en- 
ergy broadening is rendered undetectable by the statistical 
noise of the data. 

2.2 GIS 

We analyzed the GIS events files processed by the standard 
REV2 screening. The only additional screening required for 
a few observations involved removing time intervals with 
anomalously large deviations from the mean light curve of 
the data. No dead time corrections were required for the GIS 
data because the count rates were less than 1 ct s _1 for each 
GIS detector for all the observations in Table |§| 

For consistency with the SIS data, the regions used to 
extract the GIS data were chosen to be of similar size to 
the corresponding SIS regions of a given sequence. Since, 
however, the GIS+XRT PSF is somewhat larger than that of 
the SIS+XRT the extraction regions for the GIS are usually 
~ 20% larger. 

The standard observations of blank fields were used to 
estimate the background for the GIS data. The principal 
motivation for using these templates for GIS data is that 
the instrumental background increases rapidly with increas- 
ing off-axis angle. We performed the rise-time filtering (i.e. 
GISCLEAn) on the standard GIS templates as required to 
match the standard screening. In Table ^ the background- 
subtracted exposures and count rates are listed for each GIS 
observation. 

The RMFs for the GIS2 and GIS3 are equivalent (and 
are assumed to be independent of time), and thus for all 
GIS spectra we use the GIS RMFs gis2v4_0.rmf and its twin 
gis3v4_0.rmf obtained from the HEASARC archive. Using 
these RMF files an ARF file was generated for each obser- 
vation. 

2.3 Comments on individual groups 

The observations of several of the groups warrant further 
discussion. 

HCG 51: The SIS were operated for ~ 45 ks in 2-ccd mode 
and for ~ 27 ks in 1-ccd mode. We separately reduced the 
data for each mode and then combined the source and back- 
ground spectra, RMFs, and ARFs using the standard soft- 
ware (i.e. with FTOOL addascaspec vl.27). 
MKW 9: We adopted slightly stricter screening for the SIS 
to help compensate for the degradation in data quality aris- 
ing from the warmer CCD temperatures: i.e. BR_EARTH 
> 40 and RBM < 50. 

NGC 533: Each SIS for both the 1994 and 1996 observa- 
tions were unfortunately quite warm for most of the expo- 
sures. The maximum temperatures listed in Table repre- 
sent a compromise between a desire to (1) maximize data 
quality, (2) maximize the acceptable exposure, and (3) to 
best match the temperatures of the RDD maps available 
near the dates of the observations. The maximum tempera- 
tures chosen significantly reduced the effective exposure. As 
a result, since the maximum temperatures still indicate that 
the data quality are not optimal, we followed the advice in 
the ASCA ABC guide and increased the acceptable number 



© 0000 RAS, MNRAS 000, 000-000 



6 D. A. Buote 



of events by raising the threshold for pixel rejection to 100 
for both SIS of the 1996 observation; for the 1994 observa- 
tion we raised the threshold to 150 for the SISO data and to 
200 for the SIS1 data. 

Approximately 6 ks of 2-ccd mode data were also taken 
during 1994. These data were not analyzed since much of 
the group flux was lost due to the unfortunate placement of 
the group center near the top edge of the default chips SISO 
chip 1 and SIS1 chip3 (an issue more serious for the SIS1 
data). Finally, there is a noticeable point-like source ~ 7' 
from the group center. As a result, we limited the sizes of 
our extraction regions, particularly the GIS because of its 
larger PSF, to reduce contamination from this source. 
NGC 2563: Approximately 7' from the group center there 
is a point-like source. 

NGC 4104 (MKW 4s): There were two observations 
taken in 1996. One of the observations only yielded ~ 3ks of 
useful exposure when restricting the maximum SIS temper- 
atures to — 60°C and thus was not included in our analysis. 
NGC 4325: The shape of the SISO spectrum for energies 
~ 3-5 keV is quite sensitive to the size of the region cho- 
sen. As the region size is increased from our chosen r = 3' 
the spectrum in that energy range flattens significantly. It is 
possible that an unresolved hard source(s) contaminates the 
spectrum at larger radii. More likely is that the SIS back- 
ground tem plate is inadequate for this 1997 observation (see 
section 2.1). To reduce the sensitivity to either of these is- 



calculated (Liedahl et. al 1995). Although there exist other 



sues we choose the minimum region size appropriate for the 
PSF. This choice is also justified because (1) larger regions 
centered on the emission peak go off the edge of the SIS, 
and (2) the emission is very centrally peaked. 
NGC 5129: There is a fairly bright point-like source ~ 7.5' 
from the group center. 

NGC 5846: We found that the standard REV2 events files 
produce an SISO image having a noticeable excess stream of 
light on the SISO chip 0. As a result, we screened the data 
ourselves and chose slightly stricter screening criteria: i.e. 
BR_EARTH > 40, 30 for the SISO, SIS1 and RBM < 50 for 
both SIS. 

NGC 6329: The GIS2 data between 1-1.5 keV dip signifi- 
cantly, an effect which is much less pronounced in the GIS3 
data and not observed with the SIS data. We were unable 
to affect this behavior in the GIS2 data by changing the 
screening criteria. Since this dip is likely to be an instru- 
mental effect we exclude the GIS2 data from analysis. 
RGH 80: There is a fairly bright source ~ 7' from the 
group center. We limit the sizes of the regions to reduce 
contamination from this source. 



3 BROAD-BAND SPECTRAL FITTING 
3.1 Preliminaries 

3.1.1 Models 

Since the X-ray emission of the groups in our sample is dom- 
inated by hot gas, we use coronal plasma models as the ba- 
sic component of our spectral models. We use the MEKAL 
plasma code which is a modification of the original MEKA 
code (Mewe, Gronenschild, & van den Oord 1985; Kaastra 
& Mewe 1993) where the Fe L shell transitions crucial to 
the X-ray emission of ellipticals and groups have been re- 



plasma codes, most notably the popular Raymond-Smith 
code (1977, and updates), we have shown previously (B99 
and BF) that the MEKAL code provides a better descrip- 
tion of the ASCA emission of the brightest elliptical galaxies 
(which are similar to groups) than does the current version 
of the Raymond-Smith code. 

The superiority of the MEKAL code is expected because 
although both the MEKAL and Raymond-Smith codes are 
identical in their treatme nt of the ionization balance as given 
by A rnaud & Raymond (1992) for Fe and Arnaud & Rothen- 
flug ( 1985[ ) for the the other elements, the Raymond-Smith 
code has many fewer lines than does MEKAL, and, more 
importantly, it does not incorporate the im proved calcula- 
tions of the Fe L shell lines by Liedahl et al. (1995). Despite 
these differences, qualitatively similar results for ellipticals 
are usually obtained from the different codes; we refer the 
reader to section 3.2.2 of B99 for a detailed comparison. For 
our present investigation we shall employ only the MEKAL 
code. 

We account for absorption by our Galaxy using 
the photo-electric absorption cross sect ions according to 
Balucihska-Church & McCammon (1992). The absorber is 



modeled as a uniform screen at zero redshift. The column 
density is generally allowed to be a free parameter to indi- 
cate any additional absorption due to, e.g., intrinsic absorb- 
ing material, calibration errors, etc. Note, however, that the 
column density obtained in this manner only estimates the 
excess absorption because the actual value depends on the 
details of the absorber; e.g., a patchy intrinsic absorber re- 
quires factors for partial covering and the source redshift. 

We take solar (p hotos pheric) abundances according to 
Anders & Grevesse (1989) which give an Fe abundance of 
4.68 x 10 -5 relative to H by number. Since the meteoritic 
solar abundances (Fe/H of 3.24 x 10 _J by number) are prob- 
ably more appropriate, the Fe abundances we quote sh ould 
be raised by a factor of 1.44 (Ishimaru & Arimoto 1997). We 



continue to use the photospheric abundances in this paper 
to facilitate comparison with previous studies. 

The simplest spectral model we consider is that of 
isothermal gas which we refer to as the "IT" model. This 
model is expressed as phabs x MEKAL, where "phabs" rep- 
resents the absorption and "MEKAL" the thermal plasma 
emission. The free parameters of this model are the column 
density, Ah, the temperature, T, the metallicity, Z, and the 
normalization (emission measure), EM. We focus on models 
where only the Fe abundance is allowed to be a free param- 
eter and the abundances of all other elements are tied to 
Fe in their solar ratios; i.e. in this case Z indicates the Fe 
abundance. Although the Fe abundance is by far the most 
important for spectral analysis of the groups in our sample, 
we also explore models that allow the relative abundances 
of other elements to be varied separately from Fe. 

For the next level of complexity we consider models 
with two temperatures (2T) which may be expressed as 
phabs c x MEKAL C + phabs h x MEKAL h . For these "two- 
phase" models the first component represents the "colder" 
gas and the second component the "hotter" gas. The temper- 
atures and columns of each component are free parameters. 
Similar to the IT case, we focus on models where the abun- 
dances of all the elements of a plasma component are tied to 
Fe in their solar ratios; i.e. only the metallicity, Z, is a free 
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parameter. However, we tie together the abundances of the 
hotter component and those of the colder component. (In no 
case did we find that relaxing this constraint improved the 
fits.) In addition we found it useful to explore a model with 
three temperatures (3T) for some galaxies which is just the 
obvious extension of the 2T case. 

We also consider a multiphase cooling flow model (CF) 
which provides a simple model for gas emitting over a range 
of temperatures ( lohnstone et al. 1992 ) . This model assumes 



gas cools continuously at constant pressure from some up- 
per temperature, Tmax- The differential emission measure is 
proportional to M/A(T), where M is the mass deposition 
rate of gas cooling out of the flow, and A(T) is the cooling 
function of the gas (in our case, the MEKAL plasma code). 
It should be emphasized that this is arguably the simplest 
model of a cooling flow with mass drop-out. The advantage 
of this particular model is that it is well studied, relatively 
easy to compute, and a good fit to several ellipticals (e.g. 
BF and B99). 

Since gas is assumed to be cooling from some upper tem- 
perature T max , the cooling flow model requires that there be 
a reservoir of gas emitting at temperature T max but is not 
participating in the cooling flow. To accommodate this sce- 
nario we define our standard multiphase cooling flow model 
to have two components, phabs c x CF + phabsh x MEKALh, 
where the upper temperature of the cooling flow, T max , and 
the temperature of the ambient gas, Th, are tied together in 
the spectral fits. We denote this model by "CF+1T"; these 
CF+1T models thus have one less free parameter than do 
the 2T models. We also consider the cases where (1) only 
the CF component is present, and (2) a plasma component 
with variable temperature is added. 

Finally, for completeness we also consider adding a high- 
temperature bremsstrahlung component (BREM) to these 
models to account for possible emission from discrete sources 
in the central galaxies. Such a component is not expected to 
be important for groups because of their larger ratios of X- 
ray to optical luminosity with respect to isolated ellipticals 
(e.g. Canizares et al. 1987). 

3.1.2 Calibration and fitting issues 

The SIS provides substantially better constraints on the 
spectral models than the GIS because the superior energy 
resolution and larger effective area near 1 keV of the SIS 
allow better measurement of the portion of the spectrum 
where the bright, temperature sensitive Fe L shell lines dom- 
inate the X-ray emission of groups with average temper- 
atures around 1-1.5 keV (e.g. see section 3.1 of B99). In 
addition, the SIS covers energies down to ~ 0.5 keV as op- 
posed to ~ 0.8 keV for the GIS. Since there exist calibration 
problems and uncertainties with both the SIS and GIS (see 
ASCA GOF WWW pages), particularly for energies below 
1 keV, we have modified the spectra and fitting procedures 
accordingly. 

As mentioned in section |2.l| radiation damage has re- 
duced the sensitivity of the SIS for energies below 1 keV 
and has become increasingly important over the past two 
years. Consequently, for observations taken before 1996 we 
set the minimum SIS energy to 0.55 keV and raise this to 
0.65 keV for later observations. We chose this date based 
on our previous analysis of the deviations in the SISO and 



SIS1 spectra of the bright elliptical NGC 4636 observed in 
December 1995 (see Figure 5 of B99) . 

The GIS has been calibrated usi ng observations of the 
Crab nebula (Fukazawa et al. 1997). There are some re- 
ported problems with the gain of the GIS, although the en- 
ergy dependence of these problems have not yet been re- 
ported (see ASCA GOF WWW pages). Hence, we initially 
set the lower energy limit for the GIS spectra t o 0.8 keV after 
inspection of the figures in Fukazawa et al. (1997). (Below 
we revise this choice.) 

We decided to analyze SIS and GIS data below 9 keV 
(similar to B99) because at the highest energies the cali- 
bration err ors associated with the XR T become increasingly 
important (Gendreau & Yaqoob 1997), and the background 



level generally dominates the signal (thus amplifying any 
errors in the background determination). 

Errors in the calibration of the energy scale of the SIS 
and GIS result in energy offsets of a few eV for the SIS and 
~ 50 eV for the GIS (see ASCA GOF WWW pages). To 
compensate for these systematic errors we follow B99 and 
allow the redshifts for all detectors of the SIS and GIS to be 
free parameters. Note that in practice first we determine the 
best fitting spectral model with the redshifts of the detec- 
tors all fixed at their redshifts in Table |l|. Then the redshifts 
are allowed to be free to determine if there is any final im- 
provement in the fits. 

All spectral fitting was p erformed with the software 
package XSPEC ( [Arnaud 1996] ). We used the exact MEKAL 
models in XSPEC (i.e. model parameter "switch" set to 0), 
not those derived from interpolating pre-computed tables. 
For all fits we used the \ 2 method implemented in its stan- 
dard form in XSPEC. 

In order for the weights to be valid for the x 2 method 
the PI bins for each source spectrum must be re-grouped so 
that each group has at least 20 counts. This is generally re- 
garded as the standard binning procedure. The background 
files also typically have at least 20 counts when their en- 
ergy bins are grouped similarly to the source spectra. The 
background-subtracted count rate in a particular group can 
be small, but the uncertainties in the source and background 
are correctly propagated by XSPEC to guarantee approxi- 
mately gaussian statistics for the statistical weights (K. Ar- 
naud 1998, private communication). To optimize S/N we 
further rebinned the spectra into PI groups having a mini- 
mum of 40 and 80 counts respectively for the SIS and GIS 
so that the energy groups do not greatly over-sample the 
energy resolution of the individual detectors. However, we 
emphasize that (like B99) we find that the derived model 
parameters and the relative quality of the fits when using 
this courser binning are nearly identical to those obtained 
when using the standard binning. 

Hence, with the SIS and GIS spectra so prepared ini- 
tially we fitted the models described in the previous section 
jointly to the available SISO, SIS1, GIS2, and GIS3 data for 
each group. By keeping the data of these detectors separate 
we were able to assess the importance of calibration differ- 
ences between the various detectors. In particular, we al- 
lowed the column densities and redshifts of each detector to 
be fitted separately to allow for differences in the low energy 
responses and the energy scales of the various detectors. 

The procedure we follow for spectral fitting is as de- 
scribed in BF and B99. That is, we begin by fitting a IT 
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model jointly to the spectra; i.e. a single MEKAL model 
modified by Galactic absorption (see Table |l|) . The free pa- 
rameters are the temperature, metallicity, and normaliza- 
tions of each detector. We then examine whether allowing 
Ah to be free significantly improves the fit. For the 2T model 
we add another temperature component to the best-fitting 
IT model with the column density initially tied to the value 
of the first component; the redshifts and abundances are 
also tied to those of the first component. Once the new best 
2T model is found we allow A^ to vary separately from 
and fit again to obtain a new best fit. A similar procedure 
is followed for the other multi-component models. 

Upon fitting IT models to the spectra of each group we 
found that the column densities for each detector are con- 
sistent in essentially every case within their la errors; i.e. 
the statistical errors on the derived Ah are larger than the 
systematic errors due to calibration problems^]. This same 
agreement between detectors is found for JVh on both com- 
ponents of 2T and CF+1T models for seven of the groups 
in our sample. 

For the other five groups we found that £ 1 for 

the SIS data but N^/N^ < 1 for the GIS data. This occurs 
because the SIS data, which dominates the fits, generally 
requires both model components to have temperatures near 
1 keV. Because the GIS is not as sensitive to energies near 
1 keV, the GIS would rather have one component near 1 
keV and the other component with a much higher temper- 
ature. To accommodate this desire the GIS raises the value 
of A"h to a very large value (~ (several) x 10 22 cm" 2 ) which 
suppresses the contribution of the hotter temperature com- 
ponent below ~ 2 keV. Then the GIS raises the emission 
measure of this hotter component to better fit the higher 
energies. (Note that the improvement in \ 2 when allowing 
the column densities of the different detectors to be varied 
separately is always substantially less than the improvement 
of, e.g., a 2T model over a IT model.) 

Thus, in order to guarantee qualitative consistency be- 
tween the different components of the multi-temperature 
models it is necessary to tie together the column densities 
between the different detectors. This requirement is further 
justified by the agreement of the column densities of the dif- 
ferent detectors for the IT models. When tying A^h between 
different detectors the measured value simply represents an 
average. Since the constraints on such an average Ah pro- 
vided by the GIS data are dwarfed by the SIS data, and 
there exist differences in the calibration of the GIS2 and 
GIS3 below 1 keV that are currently not understood (see 
Figure 4 of Fukazawa et al. 1997), we raise the minimum 
energy of the GIS and exclude energies below 1 keV as also 
done in B99. (We do not notice any significant effect as a 
result of this choice.) 

We also find no significant differences in the measured 
redshifts between the SISO and SIS1 or the GIS2 and GIS3 
for any of the groups; i.e. only the redshift differences be- 
tween the SIS and GIS are important. Hence, it is sensible 



t Since the differences in TVh between the SISO and SIS1 are 
insignificant the less accurate RDD correction applied to the 
BRIGHT mode data of some of the observations is unimportant 
and thus is outweighed by the gain in exposure time (see section 
2.1). 



to further reduce the number of free parameters and tie the 
together the redshifts of the SISO and SIS1 and similarly the 
GIS2 and GIS3. 

Since we intend to tie together Nn between the dif- 
ferent detectors as well as the SIS0/SIS1 and GIS2/GIS3 
redshifts, we have eliminated most of the value in analyz- 
ing separately the individual detectors of the SIS and GIS. 
Since S/N is a premium for our analysis we follow B99 and 
instead analyze the summed SIS=SIS0+SIS1 and summed 
GIS=GIS2+GIS3 data. The small loss of information result- 
ing from averaging the SISO and SIS1 response matrices and 
the GIS2 and GIS3 ARF files is outweighed by the gain in 
S/N and the increase in the speed of the spectral fitting 
afforded by the reduction in the number of data sets. 

Since combining the data sets increases the number of 
counts per PI bin, we slightly increase the PI binning criteria 
to a minimum of 50 counts for the summed SIS data and 
to 100 counts for the summed GIS data. Note that similar 
to B99 we find that the best-fitting parameters obtained 
from fitting the summed SIS and summed GIS data differ 
insignificantly from those obtained when fitting all detectors 
separately (with columns tied between the detectors). 

Finally, we emphasize that obtaining the global min- 
imum for the multicomponent models often requires some 
effort (see BF and B99). Often the fits land in a local mini- 
mum upon first adding another component to a model. The 
typical case is illustrated as follows. Starting with the best 
IT model the 2T model is constructed by adding a second 
temperature component with the abundance of the second 
component tied to the first component and 7h set to some 
value greater than 1 keV. Initially one finds that the best- 
fitting abundance is small and similar to the IT case while 
Th <; 3 keV. By resetting the abundance to a large value 
(e.g. 2 solar) and Th to a smaller value near 1 keV a new 
deeper minimum is generally obtained. 

For example, the initial 2T fit to NGC 5129 gives a value 
of x 2 = 41.3 (32 dof) which is an improvement over the IT 
value of 47.8 (36 dof). The best fitting abundance is OAAZq 
which is similar to the IT value of 0.13Zq, and the best- 
fitting Tk is 100 keV which is the maximum allowed by the 
MEKAL model. However, upon resetting Z and Th following 
the above prescription the fit improves further to \ 2 = 36.1 
and gives best-fitting values Z = 0.69Zq and Th = 1.46 
keV wh ich are significantly different from the initial fit (see 
section 3.2). 



Therefore, after a fit is completed we always reset a 
subset of the free parameters (especially Z and Th) and fit 
again several times until we are satisfied that the minimum is 
stable. Stepping through the parameters when determining 
confidence limits also was useful in assessing the stability of 
the minimum. This is not the most rigorous method to find 
the global minimum, but it is at present the most convenient 
way to do it in XSPEC. 

3.2 Results 

In sum: for each galaxy group we analyze the SIS data which 
are a sum of the SISO and SIS1 data, and similarly the GIS 
data which are a sum of the GIS2 and GIS3 data. The SIS 
data are restricted to 0.55-9 keV for observations before 1996 
and to 0.65-9 keV for observations since 1996. The energy 
range for the GIS data is 1-9 keV for all observations. Finally, 
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Table 4. Spectral Fits 



Group 

HCG 51: 
IT 
2T 
CF+1T 

HCG 62: 
IT 
2T 
OF 

MKW 9: 
IT 
2T 
CF+1T 

NGC 533: 
IT 
2T 
CF+1T 

NGC 1132: 
IT 
2T 
CF+1T 

NGC 2563: 
IT 
2T 
CF+1T 

NGC 4104: 
IT 
2T 
CF+1T 

NGC 4325: 
IT 
2T 
CF+1T 

NGC 5129: 
IT 
2T 
CF 

NGC 5846: 
IT 
2T 
CF 

NGC 6329: 
IT 
2T 
CF 

RGH 80: 
IT 
2T 
CF 



(10 2 



0.69 
3.17 
2.92 

0.00 
1.32 
1.78 

1.38 
4.72 
2.28 

0.00 
2.94 
2.77 

0.02 
3.87 
3.40 

0.85 
4.36 
3.37 

0.12 
2.73 
2.23 

0.96 
3.17 
3.42 

o.oo: 

2.62] 
2.29] 

1.79 
2.69 
1 



0.22 
0.20 
+ 1.95 
3.17 
+0.94 
2.09 



+0.27 
0.00 
1.3<i 
1.32 
+0.41 
0.49 



+0.44 
0.41 
1.98 
2.97 
1.04 
2.28 



+0.60 
0.00 
+2.16 
2.94 
+0.84 
2.77 



+0.45 
O.02 
+2.07 
3.10 
+2.21 
3.40 



+0.47 
0.41 
+2.03 
4.36 
+ 0.99 
2.19 



+0.29 
0.12 
+2.37 
2.73 
+0.53 
0.54 



+0.52 
0.43 

0.80 

0.86 
+0.52 
0.67 



11.00 
1.48 
2.61 
1.48 



+0.74 
0.58 
+0.57 
0.64 
+0.64 
-0.37 



-,+0.28 
-0.00 
-1.45 
-0.42 
i 79 +0.69 



o.oo: 

0.42 4 



0.19 
2.42 
2.03 



■0.1!) 

-1.04 
-1.22 
+0.44 
0.44 



A, 



2 ) (10 2 



h 
H 

cm" 



0.77 
0.59 



+0.51 
0.41 
+0.99 

0.59 



1.44" 1 



3.14" 1 



, + 1.20 
-1.23 

2.81±i;£ 



9 09 + 1.34 
-2.04 
,+3.42 
-0.73 



0.73"* 



0.00 
0.00 



+0.63 
0.00 
+1.40 
0.00 



0.40 

0.00 



+2.18 

0.40 
+4.97 
0.00 



1.12 
0.00 



+ 1.17 

0.98 
+3.17 
0.00 



0.22 
1.42 



+0.92 
0.22 
+ 1.08 
1.07 



0.00 



+4.21 



0.00 



+0.79 



tied 



0.00 



+0.42 
0.00 



(keV) 



1.35 
0.68 
1.54 

1.07 
0.69 
1 



+0.05 
0.05 
+0.15 
0.13 
+0.17 
0.12 



+0.04 
0.04 
+0.06 
0.05 
+0.18 
-0.12 



,+0.15 
-0.14 
-+0.15 
-0.12 
9 49+0.59 
-0.33 



2.16: 

0.65" 1 



1.26 
0.63 
1.56 



+0.06 
0.09 

+0.09 
0.10 

+0.19 
0.31 



-1 ic+0.04 
1,lo -0.05 
+ 0.19 

-0.18 
+0.14 
0.09 



0.71 
1.22 



1.35 
0.80 
1.81 



+0.11 

0.10 

+0.18 

0.13 
+0.52 
0.30 



9 ic+0.13 
z ' lo -0.13 
,+0.23 
-0.12 
-1+0.32 
-0.40 



0.76: 
3.10"' 



0.85 
0.66 
0.97 

0.92 
0.68 
1.37 



+0.04 
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Best-fitting values and 90% confidence limits on one interesting parameter (Ax 2 = 2.71) are listed for isothermal (IT), 
two-temperature (2T), and multiphase cooling flow models (CF or CF+1T) jointly fit to the SIS and GIS data; the values of \ 2 i the 
number of degrees of freedo m (dof). and the y 2 null h ypothesis probability (P) are also given. Abundances are quoted in units of their 
(photospheric) solar values (Anders & Grevesse 1989); i.e. Fe abundance is 4.68 X 10 relative to H. For models with multiple 



temperature components the abundances of the hotter component are tied to the corresponding abundances of the colder component 
during the fits. For CF+1T models is tied to T c during the fits. The emission measures (EM) are quoted in units of 
10~ 17 n e n p V/47rD 2 similar to what is done in XSPEC except for the CF compo nents where instead EM C refers to M in solar masses 
per year using the distances described in the notes of Table |lj See section 3.1.1 for further explanation of the models. 
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Table 5. Three-Component Spectral Fits for NGC 5846 

AT* N* Tj T 2 T 3 Z 

(10 21 cm- 2 ) (10 21 cm' 2 ) (keV) (kcV) (keV) (Z ) EMi EM 2 EM 3 

3T: (x 2 = 81.0, 72 dof, P = 0.22) 

n -0+2.42 . 4 n+1.30 nm +0.09 n7fi+ ai ° o 97 +2.86 o oc,+oo n 14 +0.52 n rrc+0.82 n1fi +0.13 

0-4»_0.48 4 - 4i -0.92 U - 5U -0.18 U - 7 °-0.08 d - i7 -0.88 d - 8d -2.53 U - 14 -0.13 U - 55 -0.48 °- lb -0.06 

CF+2T: (x 2 = 82.7, 73 dof, P = 0.20) 
2.23±°:~ B.72±|;»3 0.77i»;JJ tied 3.72±^ 9g 2.00_ .92 5.70±J;«1 0.66±°;^ 0.20±g;iJ 

Best-fitting values and 90% confidence limits on one interesting parameter (Ax 2 = 2.71) are listed for three-component models jointly 
fit to the SIS and GIS data of NGC 5846; see Table Q for additional explanation. The column density of the third component is tied to 
Nj, for the 3T model and is fixed to the Galactic value for the CF+2T model. For the cooling flow model (1) T2 is tied to the upper 
temperature of the cooling flow (Ti) and (2) EMi is M. 



the PI bins are regrouped so that each group has a minimum 
of 50 counts for the SIS and 100 counts for the GIS. 

Table | lists the results of fitting IT, 2T, and cool- 
ing flow (CF and CF+1T) models jointly to the SIS and 
GIS data of eac h grou p. Comments on individual groups are 
given in section 3.2.E . 



3.2.1 Isothermal models (IT) 

In Figure |l| we plot the best-fitting IT models and residuals 
corresponding to the results in Table ^| for six of the groups 
possessing among the best S /N data in our sample. Only the 
SIS data is shown because of space considerations and since 
the GIS residuals are not nearly so prominent to the eye 
owing to the poorer energy resolution of the GIS. (We refer 
the interested reader to B99 for similar plots of the bright 
elliptical NGC 1399 which show GIS data and residuals.) 

In terms of the \ null hypothesis probability, P, we find 
that for most of the groups the IT models yield fits that are 
of poor quality or marginal quality (P Si 0.1). Inspection of 
Figure hi shows that the poor fits are largely the result of 
the inability to fit the SIS data near 1 keV. These groups all 
display similar residuals near 1 keV which are analogous to 
those found for the brightest elliptical galaxies by BF and 
B99. Notice that NGC 2563 has prominent 1 keV residuals 
even though P = 0.22. 

Some of the groups have significant residuals at higher 
energies such as HCG 51, NGC 4325, and NGC 5846. The 
residual pattern above ~ 2 keV for NGC 5846 clearly indi- 
cates e missio n from a higher temperature component (see 
section 3.2.2). However, for the other systems the residu- 



als at higher energies probably indic ate errors in the back- 
ground subtraction (see section 3.2.5). We mention that such 



remaining errors in the background are not very important 
for our current models because the results do not change 
appreciably if the energies above ~ 3 keV are excluded from 
the fits; i.e. the SIS data near 1 keV dominate the fits. 

Our conclusion that the IT models are generally poor 
fits to the ASCA spectra of these groups is qualitatively 
consistent with the \ 2 values (an d thei r implied P-values) 
listed in Table 2 of Davis et al. dl999| ) for several of the 
groups. (The poor quality of the IT fits is not discussed 
by Davis et al.) For other groups, especially HCG 62, NGC 
4104, and RGH 80, we obtain substantially smaller P values. 
The ability of our spectra to better discriminate between 
models is not surprising because our spectra are derived from 



much smaller spatial apertures where the S/N is optimal 
and contamination from unresolved sources is significantly 
less important. 

Most of the fitted column densities are consistent with 
the Galactic values (Table ^) when considering the SIS cal- 
ibration uncertainties and, more importantly, the statistical 
errors. Only MKW 9 and especially NGC 5846 show sig- 
nificant excess absorption, although interpretation in these 
cases is difficult because the fits are of poor quality. In con- 
trast, some galaxies have best-fitting column densities be- 
low the Galactic value which are only marginally consis- 
tent within the estimated 90% errors with the Galactic val- 
ues (e.g. HCG 62, NGC 1132, NGC 6329). For example, 
in Figure |l| the IT model for energies less than 1 keV lies 
clearly below the data for HCG 62 even with the best-fitting 
JVh = 0. Considering that the SIS calibration errors should 
translate to a measured excess Ah such examples give addi- 
tional evidence against the suitability of the IT model. 

The metallicity derived for most of the groups is very 
sub-solar: for the entire sample we obtain a mean and stan- 
dard deviation (Z) = 0.29 ±0.12 Zq. These low abundances 
are consiste nt wit h the previous study of these groups by 
Davis et al. ( [l99§ ) who fitted IT models to the ASCA spec- 
tra. Similar low Fe abundances are obtained by Fu kazawa et 
al. fll996] ) for HGC 51, by Finoguenov etal. ( |l999| ) for NGC 
5846, and by Finoguenov & Ponman fll99S| ) for HCG 51 and 
HCG 62 within our apertures. Fukazawa et al. (1998) found 
similar low Fe abundances for HCG 51, HCG 62, MKW 
9, and NGC 4104 (MKW 4s), although their analysis ex- 
cluded most of the emission at small radii which dominates 
our apertures. Such sub-solar abundances are also consistent 
with IT models fitted to elliptical galaxies (e.g. Matsumoto 
et al. 1997 ; BF; B99). However, as we show below in sec- 
3.2.3, and Appendix |a|, in most cases these low 



tions 3.2.2 



abundances result from a bias inherent in IT models. 

We have examined models where the a-process elements 
O, Ne, Mg, Si, and S are allowed to be varied separately from 
Fe. Generally the largest improvement in the fits arises from 
letting the O and Ne abundances go to zero while lesser 
improvements result when the Si and S abundances take 
values ~ 1.5-2 Zp e . For most of the groups the improvements 
in the quality of the fits are insignificant, and for those cases 
with noticeable increase in P the improvement is far less 
than when adding another temperature component as we 
discuss below. 

Similar Si and S enhancem ents are reported by 
Fukazawa et al. ( 1998 ), Davis et al. ( h.999|) , and Finoguenov 
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Figure 1. Best-fitting isothermal (IT) models fitted jointly to the SIS and GIS data for six groups with among the best S/N data in 
our sample. These models correspond to the IT results listed in Table W. Only the SIS data are shown. 



Ponman (1999). Since (1) the a/Fe ratios are generally 



poorly constrained and are consistent with solar within the 
90% uncertainties, (2) the a/Fe ratios, like the Fe abun- 
dance, turn out to be highly model dependent (see below 
and Appendix [a]) , and (3) the IT models with variable a/Fe 
ratios are in any event still inadequate fits to the data, we 
do not present detailed results of these models. We refer 
the interested reader to B99 for detailed exploration of such 
models for four bright ellipticals for which results entirely 
analogous to the groups in our sample are obtained. 



3.2.2 Two-temperature models (2T) 

The 2T models provide significantly better fits than the IT 
models to the integrated ASCA spectra of all of the groups 
in our sample with the exception of NGC 1132. The largest 
improvements occur for the groups HCG 62, NGC 5846, and 
RGH 80 which have among the lowest average temperatures 
and the best S/N data in our sample. In Figure ^ we display 
the best-fitting 2T models and residuals for the same groups 
appearing in Figure |l[ Inspection of Figure ^ reveals that 
with one exception (NGC 5846) the addition of the second 
temperature component eliminates most of the SIS residuals 
near 1 keV in a manner analogous to the brightest elliptical 
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Figure 2. As Figure |l] but for two-temperature (2T) models. 



galaxies (BF and B99). For NGC 5846 the second temper- 
ature component largely removes the resid uals for energies 
above ~ 2 keV. We explore in section 3.2.4 whether the ad- 



dition of a third temperature component can also remove 
the other SIS residuals for NGC 5846. 

Similar to the brightest ellipticals the 2T models also 
yield values of T n ^ 2 keV which indicate that the sec- 
ond temperature component arises from hot plasma and not 
the integrated emission from discrete sources in the central 
galaxy^]. For the systems with the lowest average tempera- 



's For completeness we mention that we have also investigated 
two-component 1T+BREM models where the temperature of the 



tures we obtain EM n /EM c ~ 1 in agreement with the bright- 
est ellipticals (see B99). This ratio increases with tempera- 
ture such that EM n /EM c ~ 4 for the hottest groups in our 
sample (MKW 9 and NGC 4104). 

The results listed in Table ^ for the 2T models apply to 
fits that allowed the column densities on each component to 



BREM component is set to a large value ~ 10 keV appropri- 
ate for a discrete component as indicated by some studies (e.g. 
Matsumoto ct al. 1997). For every group where the 2T model out- 
performed the IT model we also found that the 2T model fitted 
substantially better than the 1T+BREM model analogously to 
the brightest ellipticals studied by B99. 
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Figure 3. As Figure hi but for multiphase cooling flow (CF or CF+1T) models. 



vary separately (except for NGC 6329 - see section 3.2.5). 
However, the individual columns are generally poorly con- 
strained for the 2T models. (Note the lower limits on 
and upper limits on for many of the groups). The im- 
provement in the 2T models over the IT models is not due 
to the additional degree of freedom associated with allow- 
ing both iVjj and to be free parameters. The largest 
improvements are obtained for NGC 4325 and NGC 5846 
where \ 2 decreases by ~ 5 when the column densities are 
untied in the 2T fits. Finally, we note that the tempera- 
tures and abundance do not change substantially when the 
columns are varied separately. 

Despite these uncertainties we present individual con- 



straints on A^h and because for most of the groups the 
best-fitting iV^ exceeds N^, the latter of which is typi- 
cally similar to the Galactic value. This systematic trend 
is equivalent to that found for the brightest elliptical galax- 
ies, though the statistical significance is much greater for 
the latter (B99). We comment on the origin of the excess 
absorption in section ^. 

The 2T models give metallicities that are substantially 
larger than the IT models: for the entire sample we ob- 
tain (Z) = 0.75 ± 0.24 Zq; i.e. the mean is a factor of 2.6 
larger than the IT value. This systematic increase in Z has 
also been reported for the brightest ellipticals by BF and 
B99 and vividly demonstrates the sensitivity of the inferred 
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Figure 4. Best-fitting three-component (3T and CF+2T) models for NGC 5846 corresponding to those listed in Table Q. Although the 
models were fitted jointly to both the SIS and GIS data, only the SIS data are shown. 



abundance (particularly Fe) to the assumed temperature 
structure of the plasma. This systematic effect, or bias, is 
increasingly more important for the systems with lowest av- 
erage temperature; e.g., the groups MKW 9 and NGC 4104 
which have the largest average temperatures also have the 
smallest fractional increases in Z between the IT and 2T 
models. We present a detailed discussion of this fitting bias 
in Appendix 

We mention that even NGC 1132, which shows no im- 
provement in P for the 2T model, has a metallicity signifi- 
cantly larger than the IT value. Moreover, the temperatures, 
abundance, and other parameters of the 2T model agree very 
well qualitatively with those of the other groups in the sam- 
ple; i.e. NGC 1132 displays the same fitting bias as the other 
groups which strongly suggests that the similarity of its 2T 
model parameters to the other groups is not accidental. 

When allowing the a-process elements to vary sepa- 
rately from Fe the fits are not improved sub stant ially similar 
to our results for the IT models (section 3.2.1). However, 
the same abundances which improve the IT fits also im- 
prove the 2T models. Generally, most improvement occurs 
when Zo — + 0, and although the best-fitting Zsi — ■> 1.5-2 Zf c 
for some groups, the Si/Fe ratio does not differ within the 
estimated la errors from 1. 

The model dependence of the a/Fe ratios is best illus- 
trated with an example. Using IT models Davis et al. (199!:) 
infer Si/Fe enhancements for some of their groups which are 
most significant for the hottest groups MKW 9 and NGC 
4104. When allowing the a-process elements to be separately 
varied from Fe we obtain a best-fitting Zs[/Zf c = 1.4 for the 
IT model, which though significantly in excess of 1 is about 
half the value obtained by Davis et al.. This IT model also 
requires zero abundance for O and Mg and does not improve 
the quality of the fit (P = 0.1) over the case where the a/Fe 
ratios are fixed at solar. The 2T model improves slightly 
(P = 0.46) when varying the a/Fe ratios, but Zsi/Z-Fe does 
not change significantly from 1. The key improvement arises 
when letting the O abundance go to zero. 

Since the 2T models do not require Si/Fe enhancements 
for any of the groups in our sample, and the 2T models are 
clearly favored over IT models, we conclude that there is 
no evidence for such enhancements. The inferred peculiar 
O and Mg abundances inferred by 2T (and IT) models sug- 



gest that further refinement in the absorption model (see ||), 
model of the temperature structure, plasma code, or calibra- 
tion is required. We mention that if we instead add another 
temperature component (i.e. 3T model) with the a/Fe ra- 
tios fixed at the solar values typically we are able to obtain 
modest improvements in the fits of comparable magnitude 
to the models with variable a /Fe ratios discussed above (or 
larger - see below in section 3.2.4). 



3.2.3 Cooling flow models (CF and CF+1T) 

From examination of Table |i| and Figure ^ one sees that 
the multiphase cooling flow models are substantial improve- 
ments over the IT models and provide fits of comparable 
quality to the 2T models for all the groups except for NGC 
5846 and RGH 80. This similarity to the performance of the 
2T models is achieved even though the cooling flow mod- 
els have one less free parameter due to the temperatures 
of t he CF and IT components being tied together (see sec- 
tion 3.1.1). In fact, because the average temperature of the 
CF component is approximately half the upper temperature 
(e.g. Buote et al. 1998) the CF+1T models should be similar 
to 2T models with T c ~ 0.5Th which applies for the majority 
of the groups in our sample. 

For HCG 62, NGC 5129, and NGC 6329 only the CF 
component is required to produce a fit comparable to the 
2T model; i.e. for these groups (and RGH 80) adding the 
additional IT component does not improve the fits signifi- 
cantly. Note that two of these groups (NGC 5129 and NGC 
6329) have the lowest S/N data in our sample. 

In Table g the CF model is listed for NGC 5846 even 
though the CF+1T model yields a better fit with P ~ 0.05. 
The primary improvement for the CF+1T model results 
from allowing the column density on the hotter component 
to take a huge value (N^ ~ 10 24 cm -2 ) which suppresses 
the IT component for energies below ~ 2 keV and thus sim- 
ulates a higher temperature component at larger energies. 
This effect is related to the SI S and GIS column density dis- 
crepancy discussed in section 3.1.2 . Note if instead we untie 
the temperatures of the CF and IT components then a re- 
sult analogous to the 2T model is obtained. More co mplex 
models for NGC 5846 are considered below in section 3.2.4. 



The mass deposition rates (M) are listed under EM C 
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for the cooling flow models in Table ^. We obtain typical M 
values of 5-10 Mq yr~ for the groups which lie in between 
those of giant ellipticals and clusters of galaxies (see Fabian 
1994). Outliers include NGC 2563 whose small M is similar 
to bright ellipticals (e.g. B99) while the large M of NGC 
4325 is similar to clusters. The contribution of the CF com- 
ponent to the total emission measure of the CF+1T models 
is similar to, though somewhat larger than, the contribution 
of the colder component in the 2T models. This is readily 
apparent from Table ^ from comparison of the values of EMh 
for the 2T and CF+1T models. 

The column densities obtained for the CF+1T models 
agree very well with those obtained for the 2T models. Sim- 
ilar to the 2T models we also find that the improvement in 
X 2 afforded by separate fitting of A^ and A^J is insignifi- 
cant in comparison to that obtained when initially adding 
the second temperature component with JVh = N^. The 
general trend that A^ > A^j where is similar to the 
Galactic value is consistent with cooling flow models of el- 
liptical galaxies (B99) and clusters of galaxies (e.g. Fabian 
et al. 1994). 

Finally, the metallicities obtained for the cooling flow 
models systematically exceed those of the IT models and are 
similar to those obtained for the 2T models: (Z) — 0.65 ± 
0.17 Zq for the entire sample. 

3.2.4 Three-component models for NGC 5846 

NGC 5846 represents the group whose ASCA spectra are fit- 
ted the poorest by the 2T and cooling flow models: the 2T 
model gives only a marginal fit while the cooling flow model 
is unacceptable (see Table ^ . Recall that the addition of the 
second temperature component of the 2T model largely re- 
moves the residuals for energies above ~ 2 keV but does not 
significantly affect the SIS residuals near 1 keV (Figure]^). 
In this section we examine whether simple extensions of the 
previous models can eliminate the remaining SIS residuals 
for NGC 5846. 

It is possible to reduce the SIS residuals to some ex- 
tent by varying the a/Fe ratios for the 2T model leading 
to a somewhat better fit (P = 0.09). We find for such a 
model that the best-fitting temperatures, column densities, 
and Fe abundances do not change when varying the a/F e 
ratios. In contrast to the other groups (see section 3.2.2) 
we find that all of the a-process elements have best-fitting 
abundances that exceed Fe, although all the ratios except 
for Ne/Fe are consistent with a/Fe of solar within the 90% 
confidence lower limits. 

We can, however, obtain an even better fit when keep- 
ing the a/Fe ratios fixed at their solar values if instead we 
add another temperature component. In Table [| we give the 
best-fitting parameters and their 90% confidence limits for 
the 3T and CF+2T models; the corresponding best-fitting 
models and residuals are plotted in Figure ^. For the 3T fits 
we found that the column densities on the first and third 
components were similar and not well constrained individu- 
ally, and so we tied them together. Similarly, for the CF+2T 
model the column density for the third temperature compo- 
nent turned out to be small and not well constrained, so we 
fixed it to the Galactic value. 

These three-component models are formally acceptable 
and largely eliminate the SIS residuals near 1 keV similar 



to the 2T and CF+1T models of the other groups in the 
sample. The best-fitting metallicities are larger than those 
obtained for the two-component models but are consistent 
within the uncertainties. Note in Table |E] that no upper limit 
is given on Z for the CF+2T model because the current 
implementation of the CF model does not allow Z > 2Zq . 

Although the column densities and are not 
well constrained, both the 3T and CF+2T fits indicate that 
A^ > A^ in contrast to the two-component models of NGC 
5846 and of every other group in the sample. However, sim- 
ilar to the two-component models the separate fitting of the 
column densities does not improve the 3T or CF+2T model 
fits very much. For example, if the column densities of each 
component of the 3T model are tied together the result- 
ing fit is of essentially the same quality, P = 0.19, although 
some of the fitted parameters change noticeably: best-fitting 
7i = 0.32 keV, T 3 = 2.5 keV, and Z = 2.7 Zq. Nevertheless, 
it is possible that the anomalous N^/N^ ratio reflects com- 
plex absorption in NGC 5846, a galaxy with a well-known 
and prominent dust lane. 

The origin of the third temperature component is un- 
clear. It may arise in part from discrete sources in the central 
galaxy, but the temperature appears to be too low for that to 
be the only explanation. An intriguing possibility is that it 
represents an extended int ragro up gas component discussed 
by Mulchaey & Zabludoff ( 1998 ) . Since the inferred values of 
T3 may be too high for gas in hydrostatic equilibrium in the 
group potential, but are yet too small for discrete sources, 
the third component may well be a combination of the two 
effects. 

However, as noted above the 2T fits for NGC 5846 can 
be improved by varying the a/Fe ratios. In fact, we have 
found that the parameters of the 3T fits depend on the as- 
sumed a/Fe ratios. For example, let us c onsider the mete- 
oritic solar abundances of Feldman ( 1992] ) for which the Fe 
abundance is a factor of 1.44 larger than the photospheric 
case. Similar to before we find that a 3T model fits better 
(X 2 = 87.7/dof=72/P = 0.10) than the 2T model, but now 
the temperature of the second component is T2 = 2 keV. 
An even better fit is obtained by replacing the third tem- 
perature component with a BREM component; i.e. \ 2 = 
84.2/dof=72/P = 0.16 for the 2T+BREM model. For each 
of these (meteoritic) 3T and 2T+BREM models the fitted 
parameters for the lower temperature components are more 
similar to those of the other groups (Table ^]) and the bright- 
est ellipticals (B99). For example, for the 2T+BREM model 
we obtain best-fitting parameters = 1.9 x 10 21 cm" 2 , 
A/g = 0, Ti = 0.62 keV, T 2 = 1.66 keV, T B rem = 3.98 keV, 
Z = 2.79%, and BMi : EM 2 : EMbkbm = 7.4 : 2.0 : 0.17. 
Hence, using meteoritic abundances the parameters of the 
2T+BREM model for NGC 5846 are analogous to those of 
the brightest elliptical galaxies (B99), and thus the third 
spectral component (i.e. BREM) is suggestive of emission 
from discrete sources in this case. 

Finally, we mention that the other groups do not show 
as much improvement as NGC 5846 for 3T and CF+2T 
models because the SIS residuals near 1 keV are already 
most ly eli minat ed by the 2T and CF+1T models (sections 
3.2.2 and 3.2.3). The marginal fits for a few of the other 
groups are the result of other residuals (e.g. for HCG 51 at 
high energies) that are probably the result of errors in the 
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background subtraction or calibration, though errors in the 
temperature model and plasma code also may contribute. 



3.2.5 Comments on individual groups 

The spectral fitting of some of the groups require additional 
comment. 

HCG 51: If we exclude the energies above 3 keV the fits 
are improved significantly; e.g. P = 0.19 for the 2T model. 
MKW 9: The SISO and SIS1 data are very inconsistent 
below 0.7 keV such that models always lie significantly below 
the SISO and above the SIS1 data. As a result we excluded 
SISO and SIS1 data below 0.7 keV. 



NGC 533: Trinchieri, Fabbiano, & Kim (|P997j) fitted a 2T 
model to their ROSAT PSPC spectra accumulated within a 
radius of 6'. Although they tied together the column den- 
sities on both components and fixed the Fe abundance at 
solar, they also found that the 2T model fitted better than 
an isothermal model and obtained values for T c and Tk in 
good agreement with ours. 

NGC 1132: We note that the best-fitting I T abu ndance of 
O.58Z obtained by Mulchaey & Zabludoff fll998| ) using the 
Raymond-Smith code is about twice the value listed in Tabic 
^. We reproduce their result when using the Raymond-Smith 
code. 

NGC 5846: Our IT and CF models are consistent with 
those published by BF, but the values of Th and Z for our 
2T model are larger than those obtained by BF. Roughly 
equal contributors to these differences are that BF (1) set 
N u = Ah and ( 2 ) did not include the CIS data in the fits. 
NGC 6329: Since A H is unconstrained within the esti- 
mated 90% uncertainties we tied together the column den- 
sities in the 2T fits. 



4 DISCUSSION 

4.1 Comparison with previous work 

The results we have obtained for the metal abundances in 
the hot gas of poor groups of galaxies differ markedly from 
those obtained by previous X-ray studies. In this section we 
discuss more generally how our findings compare with pre- 
vious studies. As a consequence of this discussion we present 
a synthesis of the properties of the hot gas in galaxy groups 
below in section 



4.5 



We have focused our analysis of poor galaxy groups on 
the integrated ASCA spectra accumulated within apertures 
of ~ 3'-5' radius (i.e. ~ 50 - 200 kpc). For the groups in 
our sample these apertures typically enclose ~ 50 per cent 
of the total X-ray emission and, perhaps more importantly, 
<; 90 per cent of the flux within the projected aperture ra- 
dius, iiap, originates from within the three-dimensional ra- 
dius r = i? ap ; i.e. essentially our analysis is that of the in- 
tegrated spectra accumulated within the three-dimensional 
radius r ~ 50 - 200 kpc of each group. 

These region sizes were chosen as a compromise be- 
tween S/N and a desire to enclose a large fraction of the 
group emission. Since the X-ray surface brightness decreases 
rapidly with increasing radius the ability to constrain an 
arbitrary spectral model diminishes accordingly. Some au- 
thors have avoided the bright central region because models 



more complex than a single isothermal component are of- 
ten required for, in particular, galaxy clusters, whereas the 
lower S /N spectra at larger radii usually require only a sin- 
gle isothermal component. We prefer to focus on the regions 
which provide the strongest model constraints even though 
their spectra are expected to be more complex than at larger 
radii (because of, e.g., cooling flows). 

Previous X-ray studies of groups of galaxies that fitted 
a single IT model (e.g. Fukazawa et al. 1996, 1998; Davis 
et al. 1999) to the integrated ASCA spectra over a large 
region in every case obtained very sub-solar Fe abundances 
and for many cases obtained Si/Fe ratios in excess of solar. 
However, we have shown that the ASCA spectra accumu- 
lated within apertures of ~ 3'-5' radius (i.e. ~ 50 - 200 
kpc) strongly favor models consisting of at least two tem- 
perature components for 11 of the 12 groups in our sample. 
Moreover, the Fe abundances determined for the multiphase 
models (2T or cooling flow) are nearly solar and the a/Fe 
ratios are consistent with the solar values. As a result for 



studies like Fukazawa et al. (1996), who analyzed the accu 



mulated ASCA spectra of HCG 62 within an aperture of 5' 
radius, we conclude that their very sub-solar Fe abundance 
and Si/Fe ratio in excess of solar are clearly an artifact of 
forcing isothermal models to fit the multitemperature ASCA 
spectra. This fitting bias is examined in detail in Appendix 
lAl (The same holds for ROSAT studies such as Ponman et 
al. 1996 and Mulchaey & Zabludoff 1998 who found very 
sub-solar metallicities within radii of ~ 200 kpc of several 
poor groups.) 

The metal abundan ces ob tained by both Fukazawa et al. 

( 1998[ ) and Davis et al. (1999) also agree quite well with the 
values we obtained using isothermal models (Table Q), and 
thus are consistent with this fitting bias. However, it is pos- 
sible that a strong metallicity gradient may in part explain 
the very sub-solar Fe abundances obtained by Fukazawa et 
al. and Davis et al. who analyzed the accumulated ASCA 
spectra within radii that are typically factors of 5-10 larger 
than o urs. T his explanation is more plausible for Fukazawa 
et al. (1998) since they include emission only within an an- 



nulus of 0.1-0.4 h r ,n Mpc for each group. For example, for 

0.15±0.03Z Q within 
an annulus 70-286 h^ kpc as opposed to our value of, e.g., 



50 

HCG 62 Fukazawa et al. obtain Z-p 



+0.94 
-0.38^0 



Zr^. for the 2T model within a radius of 55 h 7 



0.99 

kpc. In order for these values to match there must be a 
steep metallicity gradient. The required gradient seems to 
be implausibly steep since (1) our multitemperature spectral 
models (2T and CF+1T) did not require different metal- 
licities on the different spectral compone nts for any group 
in our sample, and (2) below in section 4.4 we show that 



the Z-Rap correlation argues against a substantial decrease 
in metallicity out to radii of at least ~ 200 kpc in galaxy 
groups. 

It is even more doubtful that abundance gradients can 
give a full accounting of our differences with Davis et al. 
because our ~ 100 kpc apertures enclose typically ~ 50 per 
cent of the emission within the apertures used by Davis et 
al.. If we again consider the case of HGC 62, Davis et al. 



+0.03 , 
-0.04^0 



ZfD within a radius of 23' which 



obtained Z = 0.21 
only slightly larger than Fukazawa et al.'s value of Zp e = 
0.15±0.03Zq obtained from a region that excludes the cen- 
ter but includes the extended emission present in the aper- 
tures of Davis et al.. If abundance gradients are to account 
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for our differences, then Davis et al. should have obtained a 
weighted average of our metallicity (0.99±g;l|^ Q ) and that 
of Fukazawa et al.. This would amount to Z ~ O.GZq since 
our region and Fukazawa et al.'s region have approximately 
the same fluxes. This value is well in excess of that obtained 
by Davis et al.. Even in the extreme case that Z = out- 
side of our aperture, we would have expected Davis et al. to 
measure Z ~ 0.5Zq. 

Therefore, even if Zf c <C IZq at radii much larger than 
our apertures, the "Fe bias" is very likely an important con- 
tributor to the very s ub-solar Fe abundances mea sured by 
Fukazawa et al. (jl99§) and especially Davis et al. ( |1999| ). 

These same arguments also apply for the Si/F e enh ance- 
ments found f or som e groups by Fukazawa et al. (1998) and 
Davis et al. ( |L999t ) ; i.e. we find that S i/Fe enhancements 
can be explained by a bias (see section A2) resulting from 
fitting an isothermal model to an intrinsically multitemper- 
ature spectrum, however it is possible that abundance gra- 
dients partially could account for the fact that Fukazawa 
et al. and Davis et al. measure Si/Fe enhancements in some 
groups whereas our multitemperature models do not require 
Si/Fe ratios different from solar. 

It should be remarked that these previous studies have 
used the photospheric solar abu ndanc es whereas, as empha- 
sized by Ishimaru & Arimoto (1997), the meteoritic solar 
abundances are more appropriate. (Recall that we have used 
the photospheric solar abundances t o faci litate comparison 
with these previous studies - section 3.1.1.) For comparison 
we have also performed the spectral fits for the groups with 
the highest S/N in our sample (i.e. those in Figu res [l | - ^ 
using the meteoritic solar abundances of Feldman (IH392f) . As 
would be expected we obtain fits of comparable quality to 
before, but now the values of the Fe abundances are raised by 
a factor of ~ 1.44. As before we also find that the 2T models 
do not generally require a Si/Fe ratio different from (now 
meteoritic) solar within the 90 per cent confidence limits. 
However, the best-fitting values of the Si/Fe abundance ra- 
tio now tend to be 10-20 per cent below solar. Hence, within 
our ~ 100 kpc apertures we find no evidence for Si/Fe en- 
hancements, and when meteoritic abundances are used if 
anything we detect weak evidence for Si/Fe decrements. 

For several groups in our sample the temperature pro- 
files measured from spatially resolved spectral analysis of 
the ROSAT PSPC (e.g. Trinchieri et al. 1997; Mulchaey 
& Zabludoff 1998) and ASCA (Finoguenov et al. 1999; 
Finoguenov & Ponman 1999) data provide corroborating ev- 
idence for multi-temperature hot gas in groups of galaxies. 
The temperature profiles are qualitatively similar for each 
group: starting from a minimum at the center the tempera- 
ture rises sharply until reaching a maximum at radii of a few 
arcminutes and then falls gently at larger radii. These quali- 
tative features are entirely analogous to those of the bright- 
est elliptical galax ies; see, e.g., the compilation by Brighenti 
& Mathews (1997). Given the similarity of the X-ray spec- 
tra of poor groups and bright elliptical galaxies it should 
be expected that the very sub-solar Fe abundances of el- 
liptical galaxies found by most previous X-ray studies also 
arise from the fitting bias associated with the assumption of 
isothermal gas (BF; B99). 

The temperature maxima occur almost precisely at the 
edges of the apertures we use to extract the ASCA spectra 
(Table H); i.e. our apertures include the regions enclosing the 



positive temperature gradients in the groups (analogously 
to our previous studies of ellipticals) . Plausible physical ex- 
planations of these rising temperature profiles are, e.g., a 
single-phase medium resting in a hierarchical gravitational 
potential, a two- pha se medium, or a multiphase cooling flow; 
below in section [4.3| we discuss these scenarios and how they 
may be distinguished with future observations. 

B99 showed that the (positive) ROSAT temperature 
gradients of the bright ellipticals NGC 1399, NGC 4472, 
and NGC 5044 are inconsistent with the isothermal models 
obtained from fitting the integrated ASCA spectra accumu- 
lated within apertures of ~ 5' radius but are very consis- 
tent with simple two-temperature and cooling flow models. 
(Note the X-ray emission of NGC 5044 is b ette r described 
as that originating from a group - see section 4.4.) Since the 
ROSAT temperature profiles and ASCA multitemperature 
models for the groups are so similar to those of the bright- 
est ellipticals, we refer the interested reader to section 5 of 
B99 for a detailed demonstration of their consistency. This 
"equivalence" of a 2T model a nd te mperature gradient was 
also shown by Trinchieri et al. ( 1997 ) for the N GC 533 group 
using ROSAT. As mentioned in section 3.2.5 our 2T model 
for NGC 533 (Table §) agrees well with that of Trinchieri et 
al.. 

For our present discussion we emphasize that the pa- 
rameters obtained from fitting two-temperature and cooling 
flow models to the integrated spectra within a single aper- 
ture are self-consistent in the sense that the parameters ac- 
curately represent the values of the radially varying parame- 
ters within that aperture. We demonstrated this consistency 
via simulated ASCA SIS observations of the brightest ellipti- 
cals in section 5.3 of B99. Specifically, for each elliptical B99 
constructed a series of IT models from the ROSAT data as 
a function of radius within the radius of the ASCA aperture. 
B99 summed these radially varying models to form a com- 
posite emission model within the ASCA aperture and then 
simulated an ASCA SIS observation with the same exposure 
as the real observation. 

B99 demonstrated that the values of T c and Th deduced 
from fitting the simulated SIS data with two-temperature 
models agreed very well with respectively the lowest and 
highest temperatures of the input radially varying ROSAT 
models. More importantly for our present discussion, the 
metallicity obtained from the 2T fit reflected an average 
value of the metallicities of the input ROSAT models. This 
same level of agreement was also found for the cooling flow 
models. This exercise shows that the metallicities obtained 
by two-temperature and cooling flow models fitted to the ac- 
cumulated ASCA spectra within radii enclosing the regions 
of positive temperature gradient reflects an average of the 
metallicities over the radii within the aperture. (The con- 
sistency of the ROSAT and ASCA models fitted to the real 
data is discussed in section 5.2 of B99.) 

As a result, the metallicities of our 2T and cooling flow 
models obtained for poor groups of galaxies are inconsis- 
tent with those obtained by F inoguenov et al. (1999) and 
Finoguenov & Ponman fll999| ) within the radii we have 
examined (Table ^). These authors performed a spatially 
resolved spectral analysis of ASCA and ROSAT data for 
the groups HCG 51, HGC 62, and NGC 5846 (also NGC 
5044) and obtained very sub-solar Fe abundances and Si/Fe 
ratio enhancements that agree very well with the previ- 
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ous studies cited earlier in this section that simply fitted 
isothermal models to the integrated ASCA spectra of these 
groups; recall that the abundances obtained from those stud- 
ies suffer from the fitting bias described in Appendix |X| 
Within the radii of our apertures the Fe abundances ob- 
tained by Finoguenov et al. and Finoguenov & Ponman are 
Z Fc < OAZq for HCG 51, Z Fc < 0AZ Q for HCG 62, and 
Z Fc < 0.3Z Q for NGC 5846: These very sub-solar Fe abun- 
dances cannot be reconciled with the metallicities we have 
obtained for these groups using multitemperature models 
(Tables § and|). 

Since our two-temperature and cooling flow models ac- 
curately reflect the average metallicity within the apertures 
enclosing the region of positive temperature gradient, and 
since the two-temperature and cooling flow models (whose 
differential emission measures have very different shapes) 
generally give similar results for the abundances, we con- 
clude that the metallicities we have derived are robust and 
thus Z > 0.5Zq (90 per cent confidence) within radii of 
~ 100 kpc when averaged over the groups in our sample. 

In contrast, we are unable to assess the reliabil ity o f 
the ASCA deconvolution ana lysis of Finoguenov et al. ( 199!: ) 
and Finoguenov & Ponman (199E) since they do not demon- 
strate the validity of their procedure using simulated data. 
There are several assumptions made in their procedure (e.g. 
no correlation between spatial and spectral properties - see 
equation 3 of the Appendix in Finoguenov et al., the regular- 
ization assumptions in section 2 of Finoguenov & Ponman, 
and the assumption that the gas is single-phase everywhere) 
which could potentially impact their results. As a substitute 
to simulation, these authors present analogous results using 
ROSAT data which (in Finoguenov & Ponman) are obtained 
by using an independent method. Although the abundances 
obtained by Finoguenov et al and Finoguenov & Ponman for 
the ROSAT data of HGC 62, NGC 5044, and NGC 5846 are 
factors of 1.5-2 times larger than those obtained using the 
ASCA data, these systematic differences are dismissed be- 
cause they are marginally consistent within the 90 per cent 
confidence limits with the ASCA results. (Note no ROSAT 
results are presented for HCG 51.) 

The systematic differences of the ROSAT and ASCA Fe 
abundances obtained by Finoguenov et al and Finoguenov 
& Ponman indicate that systematic errors are significant 
either for their ASCA results, their ROSAT results, or both. 
There is ample reason to question the reliability of their 
method to deconvolve ASCA data because the results which 
have been obtained for groups and galaxy clusters are very 
sensitive to the deconvolution method used (e.g. Markevitch 
et al. 1998; Kikuchi et al. 1999; for a review of this subject 
see Irwin, Bregman, and Evrard 1999). Moreover, the PSF 
should not be very important for the ROSAT analysis but it 
is crucial for the ASCA analysis; recall that the radii under 
consideration are 3'-5', and the half-power diameter of the 
ASCA XRT is ~ 3', so the annuli analyzed by Finoguenov 
et al and Finoguenov & Ponman within our apertures are 
considerably smaller than the ASCA PSF. 

Hence, because Finoguenov et al and Finoguenov & 
Ponman do not demonstrate the reliability of their ASCA 
deconvolution via simulation, and they obtain systematic 
differences in metallicities between the ROSAT and ASCA 
data (with the ROSAT metallicities being systematically 
larger and in better agreement with our single-aperture anal- 



ysis) we suspect that the abundances determined by these 
authors within our apertures using ASCA data are incon- 
sistent with our results owing to systematic errors in their 
deconvolution procedure. 



4.2 Properties of the hot gas in groups of galaxies 

Therefore, very similar to the brightest ellipticals (B99 and 
BF) the X-ray spectral properties in the central regions (~ 
100 kpc) of groups of galaxies can be summarized as follows: 

(i) The hot gas within the ASCA apertures (~ 100 kpc) 
consists of at least two temperature components. 

(ii) The average temperature is maximum near the edge 
of the ASCA apertures and declines monotonically to a min- 
imum value at the center as indicated by both ROSAT and 
ASCA studies. 

(iii) The Fe abundances are ~ 0.7Zq and the cv/Fe ra- 
tios are consistent with their (photospheric) solar values. 
These abundances represent values averaged over the ASCA 
apertures and thus may reflect abundance gradients within. 
Though abundance gradients within the ASCA apertures 
are not required by our multicomponen t sp ectral models^], 

is suggestive of 



the Z-R correlation discussed in section 4.4 
Z declining with increasing radius. 

(iv) The two-component ASCA models suggest absorp- 
tion in excess of the Galactic value on the colder temper- 
ature component while the column density on the hotter 
component is generally consistent with Galactic. Hence, the 
excess absorption is concentrated at the center of the group. 

Alth ough t he m easure d excess absorption for the groups 
(sections 3.2.2] and 3.2.3) is not as significant as for the 
brightest ellipticals (B99) or for clusters of galaxies (e.g. 
Fabian et al. 1994), the impression is generally the same; 
i.e. iVn > A^h and JVjJ is similar to the Galactic value. The 
large values of jVjj for ellipticals and clusters imply large 
amounts of intrinsic absorbing material that generally ex- 
ceed the amount of cold gas inferred from HI or CO obser- 
vations of ellipticals (e.g. Bregman, Hogg, & Roberts 1992) 
and clusters (e.g. O'Dea et al. 1994). Although excess ab- 
sorption is consistent with the standard multiphase cooling 
flow models with mass drop-out (e.g. Fabian 1994), this dis- 
crepancy between X-ray observations and observations at 
other wavelengths poses a serious challenge to our under- 
standing of the hot gas in these systems. 

Finally, we note that our analysis does not probe the 
outer regions (i.e. r <; 200 kpc) of these groups, and thus it 
is possible that the hot gas in the outer regions of groups 
has very sub-solar Fe abundances and Si/Fe ratios above 
(phot osphe ric) solar as suggested in part icula r by Fukazawa 
et al. (1998) and Finoguenov & Ponman (199!:). However, we 
regard the ir e vidence as tentative for the reasons discussed 
i.e. isothermal modeling by Fukazawa et al. 



in section 



4.1 



and systematics in the ASCA deconvolution procedure by 
Finoguenov & Ponman. 



^ That is, for each group the abundances of each spectral com- 
ponent are generally consistent with each other. 
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4.3 Two-temperature vs. cooling flow models 

Similar to the brightest elliptical galaxies the ASCA spec- 
tra of galaxy groups accumulated within radii of ~ 3'-5' do 
not distinguish clearly between two-temperature and sim- 
ple multiphase cooling flow models; i.e. whether the hot gas 
emits at only two temperatures or from a range of temper- 
atures cannot be determined presently. This inability of the 
ASCA data to discriminate between 2T and cooling flow 
models is expected even if the cooling flow models are cor- 
rect, because the 2T models have an additional free param- 
eter and can approximate very well a cooling flow spec- 
trum throughout the ASCA bandpass (Buote, Canizares, 
& Fabian 1999). In this section we discuss possible physi- 
cal scenarios for the two-temperature models and at the end 
briefly mention the prospects for distinguishing 2T models 
from cooling flows with future X-ray observations. 

The most widely explored two-phase model emphasizes 
the role of a "hierarchical potential" structure in the hot 
gas of ellipticals, groups, and clusters. This model ascribes 
the central enhancements in X-ray surface brightness and 
central depressions in X-ray temperature to the fact that 
the central gas sits in the shallower potential well associated 
with the halo of the cD galaxy while the more extended hot- 
ter gas is associated with the gravitational potential of the 
surrounding group o r cluster. In part icular, for th e clusters 



of ga laxies Hydra- A (Ikebe et al. 1997) and A 1795 (Xuetal 



199S ) it has been shown that the two-phase models fit the 



ASCA data as well as multiphase cooling flow models. (Typ- 
ically cooling flow models of clusters also incorporate the in- 
fluence of the central galaxy on the gravitational potential; 
e.g. Thomas, Fabian, & Nulsen 1987; White & Fabian 1995) 
This "hierarchical potential" model has also been in- 
voked to explain the ASCA data of the Fornax clust er which 
houses the bright elliptical NGC 1399 at its center ( Ikebe ct 



al. 1996). However, this two-phase model cannot explain the 



multi-temperature structure implied by the 2T models ob- 
tained by BF and B99 for the ASCA data within a radius 
of ~ 5'. Firstly, Ikebe et al. assume both components are 
isothermal. Secondly, the extended component correspond- 
ing to the group halo contributes only a small faction of the 
emission within a ~ 5' radius whereas the colder and hotter 
temperature components obtained by the 2T models of BF 
and B99 essentially contribute equally to the emissivity. This 
also applies t o the groups in our sample for which Mulchaey 
& Zabludoff (1998) have found similar two-component struc- 
ture in the radial distribution of the X-ray emission^. 

A two-phase model which can explain the temperature 
structures implied by the 2T fits to the ASCA spectra of the 
brightest ellipticals and groups is described by B righen ti & 
Mathews (1998; 1999a). Brighenti & Mathews ( |l999a| ) as- 
sume that primordial gas at large radii is shock-heated to 
the virial temperature of the surrounding group potential 
as the gas flows in to the center during secondary infall. 
This hotter gas phase is similar to the extended component 
invoked in the "hierarchical potential" models above. In ad- 
dition, Brighenti & Mathews account for the mass ejected 



II It is possible that the third temperature component for the 
3T and CF+2T models of NGC 58 46 ma y arise in part from an 
extended group component (section p.2.4|). 



from stars in the central galaxy during the course of normal 
stellar evolution. This process continually ejects gas prefer- 
entially at the center of the system with a temperature ap- 
propriate to the gravitational potential of the central galaxy; 
i.e. the ejected gas from stellar mass-loss is cooler than the 
ambient shock-heated pri mordial gas. 

Brighenti & Mathews (1999a) have shown that this two- 
phase model can accurately reproduce the ROSAT temper- 
ature and surface brightness profiles of the bright elliptical 
NGC 4472. Since we demonstrated in B99 that the ROSAT 
temperature gradients for four of the brightest ellipticals (in- 
cluding NGC 4472) are very consistent with the 2T models 
obtained from fitting the ASCA spectra, we conclude that 
the scenario described by Brighenti & Mathews can also ex- 
plain the 2T models obtained from fitting the ASCA spectra 
of ellipticals and groups (BF; B99; this paper). 

The model of Brighenti & Mathews is a bona fide two- 
phase system whereas the above "hierarchical potential" 
model is truly single-phase (but not necessarily isothermal) 
in the sense that a single value of the density and tempera- 
ture exists at each radius. Although it has no effect on their 
calculations, it is assumed by Brighenti & Mathews that the 
stellar ejecta and ambient hot gas mix instantaneously so 
that a single-phase medium results. However, near the cen- 
ter where the injection of stellar material is most important 
this may not be the case, and the two phases may co-exist. 

In contrast, a continuum of phases is predicted to ex- 
ist in the cores of ellipticals, groups, and clusters in the 
standard multiphase cooling flow scenario (Fabian, Nulsen, 
& Canizares 1984; Nulsen 1986, 1998; Thomas et al. 1987; 
White & Sarazin 1987; see Fabian 1994 for a review). This 
scenario assumes that in the regions of highest density where 
the cooling time of the ambient gas is sufficiently short (usu- 
ally less than the assumed cluster age) gas cools rapidly be- 
low X-ray temperatures and drops out of the flow typically 
with a mass deposition profile M(< r) ~ r. The fate of the 
cooling gas remains the most serious problem for this sce- 
nario (e.g. Fabian 1994), though the cooling gas could pro- 
vide a natural explanation of the excess absorption obtained 
from fitting multi-temperature spectral models to the ASCA 
spectra of ellipticals, groups, and clusters. Although excess 
absorption is not an obvious prediction of the other mod- 
els described above, Pellegrini & Ciotti (199c) argue that 



single-phase cooling flows with partial winds can produce 
cold gas that would give rise to excess absorption. 

It is probably inevitable that aspects of both the two- 
phase models and multiphase cooling flow scenario operate 
to some extent in ellipticals, groups, and clusters: i.e. gas 
cooling out of the hot phase, cool gas ejected from stars in 
the cD, and hotter gas due to primordial gas shock-heated to 
the virial temperature of the surrounding group or cluster. 
Data from future X-ray missions will be able to distinguish 
between 2T and cooling flow m odels from spectra l fitting 
within a large spatial aperture (Buote et al. 1999). More- 
over, the vastly improved spatial resolution of the Chandra 
and XMM satellites over previous missions will allow the 
location of excess absorbing material and the spatial distri- 
bution of distinct gas phases (if present) to be mapped with 
unprecedented accuracy. 
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Figure 5. The metallicities and the emission-measure-weighted 
temperatures for ellipticals and groups obtained from two- 
temperature spectral models. The filled circles represent the 
groups from this paper while the open squares refer to the bright 
ellipticals NGC 1399, NGC 4472, NGC 4636, and NGC 5044 ana- 
lyzed by B99. The region denoted "Clusters" ind icates the typical 
area populated by rich clusters. See section 4.4 for details. 



4.4 Implications of nearly solar abundances 



In Figures 5 and 6 of Renzini ( 1997 ) it is shown that the Fe 
abundances of groups and clusters determined from X-ray 
observations are much more strongly correlated with tem- 
perature than with the luminosity Lb- For systems above 
T ~ 2 keV the Fe abundance is approximately constant at 
~ 0.3Zq. Below T ~ 2 keV there is considerable scatter, 
but most of the systems have Zp e ^ 0.2Zq . The Fe abun- 
dances of groups quoted by Renzini (1997) were obtained 
from published isothermal spectral models of the ROSAT 
or ASCA data in the literature and therefore suffer from the 
fitting bias discussed in Appendix [X| 

We can achieve a more consistent comparison using the 
IT models obtained for the groups in our sample. The aver- 
age metallicity and la range fo r the g roups in our sample de- 
duced from IT models (section ETgT ), (Z) = 0.29T0.12 Zq, 



is very consistent with those o btaine d for galaxy clusters. 
Similar to Figure 6 of Renzini (1997), Z tends to decrease 



with decreasing temperature for our IT models. This trend 
evidently continues below ~ 1 keV since, e.g., BF obtained 
(Z) = 0.19 ± 0.12 Z Q from IT fits of a sample of 20 bright 
ellipticals. 

Many have speculated that since Z correlates with T 
below ~ 2 keV, whereas Z is approximately constant at 
higher temperatures, there is likely a systematic error af- 
fecting the measurement of the Fe abundances. Probably 
the most frequent suggestion (e.g. Arimoto et al. 1997; Ren- 
zini 1997) is that the systematic trend is due to errors in the 
Fe L lines of the plasma codes, but previously we have shown 
that existing Fe L errors are not serious for temperature and 



abundance determination of normal ellipticals which should 
be the most affected (BF; B99). Our previous analyses of 
ellipticals and our current analysis of groups instead demon- 
strate that the key systematic error in the determination of 
Z results from assuming the gas to be isothermal in these 
systems. 

In Figure |E] we plot Z versus emission-measure- weighted 
temperature obtained from the 2T models for the groups 
in our sample. For comparison we also show in Figure ^ 
the region containing m ost of the cluster Z measurements 
as compiled by Renzini (1997). We note that these cluster 



metallicities were obtained from GINGA and EXOSAT data 
and thus include emission over a large portion of the clusters 
whereas our results for groups are obtained within apertures 
of ~ 100 kpc which typically contain ~ 50 per cent of the 
total X-ray emission. 

The trend is now quite different from Renzini's Figure 6 
since below ~ 2 keV we find that Z increases to ~ 0.7Zq for 
groups. Unlike the IT models there is no obvious correlation 
of Z with (T) for the 2T models; i.e. removing the "Fe bias" 
with the 2T models also largely removes the dependence of 
Z on temperature for groups. 

The 2T models of ellipticals (BF) suggest that Z may 
continue to increase for lower temperatures, though the scat- 
ter is considerable, (Z) = 0.9 ±0.7 Zq. Let us instead focus 
on the ellipticals NGC 1399, NGC 4472, NGC 4636, and 
NGC 5044 studied by B99 since they have the best avail- 
able constraints on the spectral models. Moreover, for NGC 
1399, NGC 4472, and NGC 4636 B99 included an emission 
component that accounts for emission from discrete sources 
as expected in these systems. The results for these galaxies 
are plotted as open squares in Figure ^. These objects over- 
lap in temperature with the groups and add extra scatter in 
Z. 

Since Z and (T) do not exhibit a strong correlation for 
the groups and ellipticals we investigate correlations of Z 
with other physical quantities. In Figure ^| we plot Z versus 
L x in the 0.5-10 keV band. There is a significant correlation 
in the sense that more luminous systems have lower Z. This 
trend for groups appears to connect smoothly with clusters 
(Z ~ 0.3Zq) at high L x . We mention that similar correla- 
tions exist if instead of L x we plot Z against the emission 
measure of the cooler component (i.e. EM C x 47rD 2 ) or M 
for the cooling flow models; i.e. the trend of decreasing Z 
with increasing L x (i.e. mass) can also be expressed as Z 
decreasing either with increasing luminosity of the cooler 
component of the 2T models or with increasing mass depo- 
sition rate for the cooling flow models. 

It is clear from Figure ^| that with respect to X-ray emis- 
sion the "elliptical" NGC 5044 (log 10 L x = 42.74) is more 
accurately described as a bright group which is also indi- 
cated by its large M obtained from a cooling flow model 
(B99; also see David et al. 1994). Also evident is that the 
elliptical NGC 4636 (log 10 L x = 41.53) has smaller Z than 
NGC 1399 and NGC 4472 which have higher L x . Of the 
systems analyzed by B99 NGC 4636 is the only case which 
required the ot/Fe ratios to deviate significantly from solar in 
order to achieve an acceptable fit, and it had the largest con- 
tribution from a high energy component presumably arising 
from discrete sources; i.e. the Fe abundance of NGC 4636 
should be regarded as provisional. 

The correlation of Z with L x (or with, e.g., M) still 
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Figure 6. As Figure ^except here we plot the metallicity versus (left) the (unabsorbed) luminosity in the 0.5-10 keV band (see Table 
[I]) and (right) the size of the aperture used for spectral analysis. 



possesses considerable scatter, and thus there may be other 
factors to consider. In particular, also in Figure ^ we plot Z 
versus R, where R is the radius of the circular aperture used 
to extract the SIS0 spectrum for each object (see Table |^). 
There is a clear correlation such that the metallicities com- 
puted within the largest apertures tend to have the smallest 
values. For R ~ 50-100 kpc the correlation for the groups 
has the least scatter; i.e. Z falls from ~ l^Q at R ~ 50 kpc 
to ~ 0.5Zq at R ~ 100 kpc. At larger R there is no evi- 
dence for further decline out to R ~ 200 kpc; i.e. unlike the 
correlation with L x it is not obvious from these data that 
the Z-R correlation joins smoothly with clusters. 

Nevertheless, this trend of decreasing Z with R is con- 
sistent with the presence of abundance gradients in these 
systems. Since R and L x are correlated to some extent be- 
cause brighter sources generally allow for larger extraction 
regions, we are unable at this time to determine if the funda- 
mental correlation of Z is with L x or R or whether they are 
equally important. From examination of Figure ^| we reach 
the tentative conclusion that Z decreases due to increases in 
both R and L x , but for the groups in our sample abundance 
gradients are most important for R ~ 50-100 kpc while the 
X-ray luminosity is the key determinant for log 10 L x 42.5. 
(For the ellipticals NGC 1399, NGC 4472, and NGC 4636 
these effects cannot be distinguished.) 

Renzini (1997) also shows that the gas fraction (/(,) is 



approximately constant with temperature for clusters but 
it plummets by several orders of magnitude for tempera- 
tures below ~ 2 keV. The determinations of are not ex- 
tremely sensitive to Zp e and thus our larger Fe abundances 
will not change the trend in Renzini's Figure 4 by very much; 
i.e. the well-known conclusion remains unchanged that the 
much smaller gas fractions of groups imply that a signifi- 
cant portion of their gas be expelled during their evolution 
in contrast to clusters. 

The very sub-solar Fe abundances inferred from most 
previous studies imply that Z in groups is similar to or less 



than that of clusters. Within the context of the "standard 
chemical model" which assumes star formation properties 
similar to the Milky Way, these low values of Z require that 
groups accrete significant amounts of primordial gas at late 
times after spending much of their evolution expelling gas 
as demanded by their low ft (e.g. Renzini 1994, 1997). In 
this case the gas dynamical and metal enrichment history 
of groups does not represent a smooth transition between 
ellipticals and clusters neither of which are expected to ac- 
crete large amounts of primordial gas at late times. However, 
since our new results indicate Z ~ 0.7 Zq in groups the late 
accretion of primordial gas is now not required, and thus the 
decrease of Z from ellipticals to groups and then to clusters 
can be understood in the context of the "standard model" as 
an increase in the proportion of primordial gas retained as 
the mass of the system increases from ellipticals to clusters. 

The a/Fe ratios obtained from our 2T and cooling flow 
models are consistent with the solar values and those found 
for the brightest ellipticals (BF; B99). This result taken with 
the nearly solar Fe abundances indicates that, like clusters, 
the properties of the hot gas in groups of galaxies are con- 
sistent with the Type la to Type II supernova ratios and the 
IMF being similar to those of the Galaxy (e.g. Renzini 1997; 
see also Fujita, Fukumoto, & Okoshi 1997 and Brighenti & 
Mathews 1999b). 



5 CONCLUSIONS 

We have analyzed the ASCA spectra of 12 of the bri ghtes t 
groups of galaxies previous ly stu died by Davis et al. ( 199£ ) 
and Mulchaey & Zabludoff (1999). Because of the limitations 



inherent in the ASCA PSF we restricted our analysis to 
the integrated spectra accumulated within a single circular 
aperture corresponding to radii of ~ 3'-5' (~ 50-200 kpc) for 
the groups in our sample; these apertures typically enclose 
~ 50 per cent of the total X-ray emission. Our principal 
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motivation for re-examining these data is to verify whether 
the very sub-solar Fe abundances obtained from previous 
studies are in reality an artifact of fitting isothermal models 
to the X-ray spectra as found for bright elliptical galaxies 
by BF and B99. 

Fitting isothermal models (IT) jointly to the ASCA 
SIS and GIS spectra gives results that are of poor or at 
best marginal quality for most of the groups; i.e. P ^ 0.1 
where P is the % 2 null hypothesis probability. The primary 
cause of the poor fits is the inability of the IT models to fit 
the SIS data near 1 keV. Similar to previous studies these 
IT models give very sub-solar metallicities: for the whole 
sample we obtain a mean and standard deviation (Z) — 
0.29 ± 0.12 Zq. Although allowing the a-process elements 
to vary separately from Fe often results in a/Fe ratios quite 
different from solar (especially for O and Ne), the quality 
of the fits is not improved significantly in most cases. These 
properties of the IT fits are very analogous to those of the 
brightest ellipticals (see B99). 

Two-temperature models (2T) provide significantly bet- 
ter fits than the IT models for 11 out of the 12 groups in our 
sample. The superior performance of the 2T models arises 
from their largely removing the residuals in the SIS data 
near 1 keV present in the fits of the IT models. For groups 
where the temperature of the hotter component, Tk, is < 2 
keV the 2T models give Th ~ 2T C where T c is the temper- 
ature of the colder component in the hot gas. The ratio of 
emission measures of the two components is EMh/EM c ~ 1 
for systems with the lowest average temperatures and rises 
with increasing temperature until EMh/EM c ~ 4 for the 
hottest groups. 

Although the column densities on the individual com- 
ponents are not as precisely constrained as those for the 
brightest ellipticals (B99), the impression is generally the 
same: > and is similar to the Galactic value. We 
mention that the improvement of the 2T models over the 
IT models is not the result of the extra degree of freedom 
afforded by separately varying and N^, and the tem- 
peratures and Fe abundances do not change substantially 
either. 

The metallicities of the 2T models are considerably 
larger than obtained for the IT models, (Z) — 0.75 ± 
0.24 Zq; i.e. the mean metallicity is a factor of 2.6 larger 
than the IT value. Note that the 2T models have the a/Fe 
ratios fixed at their solar values. Even if the a/Fe ratios are 
allowed to be free parameters, the 2T models models do not 
indicate Si/Fe ratios in excess of solar in contrast to previous 
studies using isothermal models (e.g. Davis et al. 1999). 

We obtain results for the multiphase cooling flow mod- 
els that are entirely analogous to those of the 2T models 
including the large metallicities, (Z) = 0.65 ±0.17 Zq. Typ- 
ically the mass deposition rates (M) are 5-10 Mq yr~ plac- 
ing the groups midway between ellipticals and clusters (e.g. 
Fabian 1994) as expected. Hence, the value of M for the 
cooling flow models, or equivalently the emission measure of 
the cooler component for the 2T models, are viable means 
of distinguishi ng g roups from ellipticals and clusters. 

we have presented an extensive compar- 



In section 4.1 



ison of our results to previous studies. We conclude that 
the very sub-solar Fe abundances and Si/Fe enhancements 
in groups of galaxies found by most previous studies (e.g. 
Fukazawa et al. 1996; Davis et al. 1999) are an artifact of 



their fitting isothermal models to intrinsically multiphase 
ASCA spectra. (We also critically evaluate some studies of 
spatially resolved spectral analysis.) These results are vir- 
tually identical to those found for the brightest elliptical 
galaxies by BF and B99. Owing to the importance of these 
results for interpreting X-ray spectra, in Appendix [X] we use 
simulated ASCA observations to explore in some detail the 
"Fe bias" and "Si bias" associated with the spectral fitting 
of ellipticals, groups, and clusters of galaxies. 

We have discussed how the properties of the hot gas in 
groups of galaxies within the central regions (r ~ 100 kpc) 
are revised in light of the 2T and cooling flow models, and 
how the properties are consistently interme diate with those 
of ellipticals and clusters (sections 4.5 and 4.4). Also con- 
sidered are physical scenarios for the 2T models and how 
they can be distinguis hed from cooling flows with future 
X-ray missions (section 4.3). Finally, we examined the im- 
plications of the nearly solar Fe abundances and solar a/Fe 
ratios for the groups and how they are consistent with pop- 
ular chemical models. 

Our results for the brightest ellipticals (BF; B99) and 
the brightest groups (this paper) demonstrate the similarity 
of the properties of the hot gas in the centers of these systems 
to cooling flow (i.e. cD) clusters (e.g. Fabian 1994). The 
different Fe abundances measured for these systems may 
reflect the sizes of the apertures used for spectral analysis; 
i.e. the physical aperture size typically used when measuring 
Zp e increases from approximately 30 kpc for ellipticals to 
100 kpc for groups up to at least 200-300 kpc for clusters. 
However, presently the role of abundance gradients cannot 
be disentangled from t hat o f the X-ray luminosity (or M) of 
the system (see section 4.4 ) . Although ellipticals and groups 
are believed to expel much of their hot gas during their 
evolution while clusters are "closed boxes" and retain all 
of their gas, these differences may not much affect the gas 
properties in their centers, and thus the strong similarity 
of the gas properties in the cores of these systems is not 
unreasonable. 

The role of abundance gradients in ellipticals, groups, 
and clusters should be clarified significantly with data from 
the upcoming Chandra and XMM satellites which combine 
high spatial and spectral resolution unlike previous missions. 
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APPENDIX A: BIASES ASSOCIATED WITH 
X-RAY SPECTRAL FITTING 

The X-ray studies of the hot gas in elliptical galaxies 
NGC 720 by Buote fe C anizares fll994p and NGC 4636 



model 3 on t he co oling flow models of MKW 9 and NGC 
4104 (section E2~3|). 



by Trinchieri et al. (1994) were the first to show that the 



ROSAT PSPC spectra of ellipticals could not distinguish be- 
tween isothermal models with very sub-solar Fe abundances 
and two-temperature models with solar Fe abun dance s. Us- 
ing simulated PSPC data Buote & Canizares (1994) fur- 



ther demonstrated that if the galaxy spectrum is intrinsi- 
cally characterized by two temperatures in the hot gas with 
solar Fe abundances but is instead fitted with an isothermal 
model then sub-solar Fe abundances necessarily result; i.e. 
the isothermal mod el g ives a strongly biased measurement 
of the Fe abundancefj. 

The strong dependence of the Fe abundance on spectral 
model has now been demonstrated using ASCA data for el- 
liptical galaxies by BF and B99 and for groups of galaxies 
in this paper. These papers also show that multiphase mod- 
els are clearly favored over isothermal models for these sys- 
tems. As a result, investigators hoping to obtain average Fe 
abundances for these objects by fitting isothermal models in 
actuality always infer strongly biased (i.e. very sub-solar) Fe 
abundances for ellipticals (e.g. Matsumoto et al. 1997) and 
for groups (e.g Davis et al. 1999). Since this "Fe bias" is 
crucial for interpreting the X-ray spectra of these systems, 
in this Appendix we attempt to elucidate this bias by using 
simulations of simple multiphase models that are consistent 
with the current ASCA observations of ellipticals, groups, 
and clusters. 

We investigate the dependence of the Fe bias on tem- 
perature by considering objects spanning the temperatures 
of the brightest ellipticals, groups, and clusters of galax- 
ies. In addition, we use ASCA simulations of these systems 
to show the dependence of the a/Fe ratios on the spectral 
model. Thus, this Appendix updates and extends the brief 
studies of t he Fe bias that appear in section 2.2 of Buote & 
Canizares ( 1994 ) and in section 3.2.2 of BF. 

The key properties of the four ref erence models we use 
for this study are listed in Table 



Al 



Each of these mod- 
els corresponds to regions of positive temperature gradient 
(e.g. see section 5 of B99). Models 1 and 2 are based on the 
two-temperature and cooling flow models of the elliptical 
galaxy NGC 1399 within r ~ 30 kpc (B99) which are also 
quite similar to the corresponding models of several of the 
lowe r-tem perature groups in this paper (see sections 3.2.2 
and 3.2.5). To explore the high-temperature groups we base 



** We note that Bauer & Bregman (1996) discuss difficulties 



with determining metallicities from PSPC observations but not 
with respect to a fitting bias resulting from ass uming an isother- 
mal model for elliptical galaxies. Kim et al. (1996) show that 



the metallicity inferred from spectral fitting of the ASCA and 
ROSAT data of the faint SO galaxy NGC 4382 is essentially zero 
for isothermal models, but it cannot be constrained (and is thus 
consistent with solar) by a two-component model with T\ ~ 0.3 
keV for the hot gas and T2 ~ 5 keV for a bremsstrahlung com- 
ponent. Although this shows the degeneracies in metallicity de- 
pending on the spectral model of a faint (low L x /Lb) elliptical, 
this does not demonstrate a bias resulting from fitting a spectrum 
with two temperatures in the hot gas (i.e. with T ~ 0.5-2 keV) 
apart from any bremsstrahlung component. 



Finally, we obtained the ASCA data of A478 from the 
HEASARC archive and have analyzed the spectra within a 
3' radius^] in a similar manner to the ellipticals and groups 
(see section [2J) : the result is listed as model 4 and is consis- 
tent with previous studies (e.g. Johnstone et al. 1992). A478 
was selected because (1) it is one of the brightest clusters 
having evidence for substantial excess absorption and (2) 
with respect to other cooling flow clusters a comparatively 
large fraction of the emission measure arises from the cool- 
ing flow component; i.e. A478 should be one of the most 
favorable clusters for examining systematics when fitting an 
isothermal model. 

We simulated ASCA observations of the models in Table 



Al using XSPEC. Because of its better energy resolution 
and sensitivity near 1 keV we focus on the SIS data. For 
consistency of presentation we use the same PI grouping 
for each simulated SIS spectrum (i.e. that of the SIS data 
of NGC 1399 - see B99). Our choices of exposure times 
for the simulations are motivated by a desire to insure (1) 
sufficiently high S/N for each case so that noise does not 
noticeably affect the fitted parameter values, (2) that the 
sizes of the residuals (in terms of \ 2 ) are similar to the 
best available ASCA data, and (3) that the simulated data 
for each model have similar S/N. Although background is 
included in the simulations, the S/N criteria stated above 
insure that its effect is only noticeable at the highest energies 
which are not of interest to our study. 



Al The Fe bias 



In Figure Al we plot the simulated ASCA SIS data of models 
1 and 2 from Table Al and also show the best-fitting isother- 
mal model and its residuals. These IT models with the a/Fe 
ratios fixed at solar give best-fitting values of Z — 0.27Zq 
for the simulated 2T model and Z = 0.32Zq for the simu- 
lated CF+1T model. (These metallicities do not differ sig- 
nificantly from the Fe abundances obtained when the a/Fe 
ratios are allowed to vary in the fits - see Table A2.) Notice 



that both the pattern of residuals and the metallicities de- 
rived from the IT fits are very similar for both models 1 and 
2 and thus the Fe bias is essentially the same for isothermal 
fits to two-temperature and cooling flow spectra. 

This Fe bias is insensitive to the presence of excess ab- 
sorption on the colder component. If instead we consider a 
model identical to model 1 but with JV£ = Njj = 5 x 10 20 
cm -2 , then the metallicity obtained from fitting the IT 
model to the simulated ASCA data has a best-fitting value 
Z = 0.25Zq which is consistent with the value obtained 
from the simulated model 1 within the la errors. Below 
1 keV the pattern of residuals for this model differs from 
model 1 such that the residual pea k below 1 keV is slightly 
less pronounced than in Figure [All 



The pattern of residuals of model 3 (Figure A2 ) is less 
complex than for the lower temperature models 1 and 2. The 
best-fitting metallicity using the IT model (a/Fe fixed at 
solar) is 0.61Zq which represents only a 40% under-estimate 

tt This radius corresponds to 342 kpc assuming a redshift of 
0.0882, H = 70 km s" 1 Mpc -1 , and U = 0.3. 
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Table Al. Models 



Model 


Type 






T c 


T h 


Z 


EM C or M 


EM h 






(10 21 cm" 2 ) 


(10 21 cm~ 2 ) 


(keV) 


(keV) 


(%) 


(see notes) 


(see notes) 


1 


2T 


2.5 


0.5 


0.75 


1.5 


1.0 


2.0 


2.0 


2 


CF+1T 


2.5 


0.5 


1.5 


tied 


1.0 


1.5 


2.0 


3 


CF+1T 


2.5 


0.5 


2.5 


tied 


1.0 


10.0 


1.5 


4 


CF+1T 


6.8 


2.2 


6.8 


tied 


0.34 


653 


22.4 



The units of emission measure (EM) are defined as in the notes to Table ^, and the mass deposition rates (M) are listed in units of Mq 
yr _1 . See text for further explanation of these models. 



Table A2. IT fits to simulated data 



odel 


7V H 


T 




Zo 








Z S 




(10 21 cm" 2 ) 


(keV) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


1 


0.74 


1.03 


0.26 


0.00 


0.00 


0.15 


0.30 


0.42 


2 


0.35 


1.16 


0.31 


0.00 


0.00 


0.27 


0.32 


0.45 


3 


0.30 


2.20 


0.66 


0.00 


1.39 


0.67 


0.71 


0.79 


4 


3.47 


5.00 


0.33 


2.34 


0.94 


0.44 


0.64 


0.55 



Best-fitting IT model parameters resulting from fits to simulated ASCA data of the models in Table Al . All elements not listed have 
abundances tied to Fe in their solar ratios. See text for discussion of the ASCA simulations. 



Two Temperature simulation: T c = 0.75 keV and T h = 1.5 keV 



Cooling Flow simulation: T - 1.5 keV 
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Figure Al. Best-fitting IT models and residuals for the simulated ASCA data of models 1 (left) and 2 (right) listed in Tabic Al 
These IT fits have the a/Fe ratios fixed at the solar values to highlight the residuals near the regions of the strong emission lines of the 
o-process elements. 
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of the true value in contrast to the ~ 70% under-estimate 
observed for models 1 and 2. The Fe bias decreases more 
dramatically as we move to still higher temperatures with 
model 4 (Figure A2). For this model we obtain a best-fitting 
value of Z = 0.30Zq using the IT model {a/Fe fixed at 
solar) which is only a 10% under-estimate of the true value. 
If we allow the a/Fe ratios to vary (Table A2) the best- 
fitting Fe abundance rises to 0.33Zq; i.e. the Fe bias for hot 
clusters is insignificant. 

A physical explanation of the Fe bias is as follows. 
For models 1 and 2 which have average temperatures ~ 1 
keV, the colder temperature components preferentially ex- 
cite emission lines in the Fe L complex below 1 keV (from 
ions Fe xvii - xxi) while the hotter temperature components 
preferentially excite the Fe L lines from ~ 1-1.4 keV (Fe xxi 
- xxiv). Since these components contribute approximately 
equally to the emission measure, the net result is that the 
spectral shape in the highly temperature-sensitive Fe L en- 
ergy region is flattened with respect to a IT model. That 
is, the spectrum of a IT model which has the average tem- 
perature of the two components of models 1 or 2 (and has 
the same Fe abundance) will preferentially excite the Fe L 
lines right at 1 keV at the expense of the other Fe L lines 
near ~ 0.7 keV and ~ 1.3 keV and thus will be more peaked 
at 1 keV than either model 1 or 2. To partially compensate 
for this, i.e. to flatten the spectral shape near 1 keV, the IT 
model decreases its Fe abundance so that the contribution 
of the flatter bremsstrahlung continuum is increased in the 
model. However, the distinctive residual pattern in Figure 



the colder temperature component. If instead we consider a 



5 x W 



Al demonstrates that the Fe L lines are sufficiently temper- 
ature sensitive (and the SIS has sufficient energy resolution) 
to thwart this attempt to successfully fit the IT model to 
the simulated two-temperature and cooling flow spectra. 

The strong temperature-sensitive Fe L emission lines are 
the key to this Fe bias. As one moves to higher temperatures 
(models 3 and 4) the Fe L lines become weaker and the 
bias decreases rapidly for average temperatures above 2 keV. 
Moreover, at higher temperatures the strong 6.5 keV Fe K 
emission line complex becomes increasingly more important 
and dominates the determination of the Fe abundance. 

A final note: since it is the variation of the spectral 
shape of the Fe L complex with temperature which gives rise 
to the Fe bias, one does not obtain a biased measurement 
of the Fe abundance using a IT model if the spectrum con- 
sists of multiple spectral components which have the same 
temperature but different Fe abundances. In this case the IT 
model yields a value for the Fe abundance that is a weighted 
average over the various components. 



A2 The Si and S bias 

Prominent residuals in the IT fits (a/Fe fixed at solar) are 
located ne ar th e st ron g emission lines of O, Ne, Mg, Si, and 
S (Figures 



Al 



and 



A2). If the abundances of these elements 



are allowed to vary separately from Fe we obtain the re- 
sults listed in Table |A^. Although in most instances the O, 
Ne, and (to a lesser extent) the Mg abundances differ sig- 
nificantly from Fe, their new values do not correspond to a 
substantial reduction of the residuals over ~ 0.5-1.5 keV. 

The IT fits for models 1 and 2 produce zero O and Ne 
abundances and Mg/Fe ratios less than solar. These results 
are very sensitive to the presence of excess absorption on 



= I.IZq, and 



model identical to model 1 but with - 
cm -2 , then we obtain Zo = 0.55Zq, 
Zyig = 0.23Zq with Zys = 0.26Zq. 

For models 3 and 4 the IT fits yield O, Ne, and Mg 
abundances very different from models 1 and 2 and from 
each other. Although the IT fit for model 3 gives zero O 
abundance its large Ne abundance is qualitatively differ- 
ent from models 1 and 2. The ratios of O/Fe, Ne/Fe, and 
Mg/Fe obtained for model 4 are all very much in excess of 
solar. Hence, the O, Ne, and Mg abundances are very sensi- 
tive to the temperature model and to the absorption model, 
and thus the systematic errors introduced when using a IT 
model to measure these abundances are more complex in 
origin than the Fe bias discussed above. 

In contrast, the Ka emission lines of Si and S are at 
sufficiently high energies so that they are not affected by 
the absorption model. The residuals near the Si (~ 1.8-2.0 
keV) and S (~ 2.4-2.7 keV) lin es a re readily apparent for 



all the models in Figures Al and A2, especially for model 4 



The IT fits with variable a/Fe ratios listed in Table A2 for 
models 1-3 indicate modest ratios of ^ 1.1 for Z^/Zyb and 
more substantial ratios of ~ 1.2-1.5 solar for Zs/Zf b . Model 
4 has the largest ratios with Zsi/Z^e ~ 2 and Zs/Zp c ~ 
1.7. Thus, the inferred values of Si/Fe tend to increase with 
temperature while the S/Fe ratio is relatively constant. 

These excess Si/Fe and S/Fe ratios inferred by IT 
models arise from two effects: (1) the colder temperature 
components present in the two-component models produce 
stronger Si and S emission, and (2) the IT model under- 
estimates the continuum at energies above ~ 2 keV. Issue 
#1 is more important for the higher temperature systems 
like clusters (model 4) while issue #2 dominates for the lower 
temperature systems like elliptical galaxies (models 1-2). 

Therefore, attempts to infer the Si/Fe and S/Fe ratios 
using isothermal models on spectra that are intrinsically 
composed of multiple temperature components will neces- 
sarily yield values for these ratios that are biased in ex- 
cess of the true values, the effect of which is most impor- 
tant for higher temperature clusters of galaxies. This bias 
may have played a rol e in t he excess Si/Fe ratios obtained 
by Mushotzky et al. (1996) and ma y in p art explain the 
trend discovered by Fukazawa et al. (1998) that the Si/Fe 
ratio rises with temperature for clusters galaxies. Although 
Fukazawa et al. attempt to eliminate any problems associ- 
ated with multiphase gas by excluding the central regions of 
their clusters, it is quite possible that their arbitrary method 
of defining the affected regions does not fully remove the ef- 
fects of multiphase gas and the resulting Si bias we have 
described. 
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